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^Studies  were  conducted  involving  the  three  phases  of 
application  of  the  AMI  and  MNDO  semiempirical  quantum  mechanical 
molecular  models  to  problems  in  chemistry:  development,  testing,  and 
use  in  specific  investigations.  As  developmental  work,  efforts  were 
made  to  obtain  parameters  for  phosphorus  valid  for  both  trivalent 
and  pentavalent  compounds .  Both  the  standard  AMI  and  MNDO 
algorithms,  and  algorithms  including  a  core  -  core  repulsion 
function  modified  to  take  into  account  the  change  in  bonding 
interactions  with  change  in  valency,  were  used  in  these  efforts. 
Contributory  to  the  testing  of  the  AMI  model,  results  for  proton 
affinities  and  deprotonation  enthalpies  were  extensively  studied. 
The  validity,  as  well  as  the  limitations,  of  using  AMI  in  studies  of 
reactions  involving  these  processes  were  evaluated.  Finally,  MNDO 
and  AMI  were  used  to  examine  two  specific  problems  of  chemical  and 
theoretical  interest,  the  condensation  reaction  of  polyketide 
biosynthesis  and  the  mass  spectral  fragmentation  of 
cis-l-nitropropene .  Results  in  the  former  study  support  the 
suggestion  that  a  crucial  factor  in  enzyme  reactions  is  the 
exclusion  of  water  from  the  reacting  system.  In  the  latter  study, 
AMI  results  were  shown  to  correlate  well  with  the  fragmentation  of 
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Chapter  1 


INTRODUCTION 


While  man  has  always  had  a  natural  curiosity  about  the 
chemical  processes  surrounding  him,  the  understanding  of  them  grew 
slowly  prior  to  the  late  18th  and  early  19th  centuries.  Dalton's 
presentation  of  his  theory  of  atoms  in  1802,  along  with  Lavoisier's 
discovery  of  mass  as  the  fundamental  quantity  in  chemical  reactions 
in  1789,  provided  the  basic  concepts  necessary  to  understand 
chemical  reactions.  Progress  in  the  determination  of  chemical 
composition  and  the  characteristic  reactivity  of  elements  and 
functional  groups,  and  development  of  bonding  theories  and  theories 
of  reaction  mechanisms  followed.  An  explosion  of  instrumentation  in 
the  20th  century  added  to  the  momentum,  providing  information 
unobtainable  from  wet  chemical  techniques.  Spectroscopic, 
crystallographic  and  diffraction  techniques  gave  new  insights  into 
molecular  structure,  bond  strength,  electron  distribution,  group 
reactivity  and  stereochemistry.  As  a  result,  a  greater  depth  of 
understanding  of  reactions  was  reached. 


i 

I 


Still,  there  remained  a  realm  of  information  unreachable  by 
either  wet  chemical  or  instrumental  techniques.  While  enthalpic  and 
entropic  properties  of  transition  states  could  be  derived  from 
experiments,  the  structures  of  the  transition  states,  and,  hence, 
the  elementary  processes  involved  in  reaching  the  critical  points  of 
the  reactions,  could  not.  Furthermore,  this  was  not  a  result  of 
insufficiently  developed  techniques,  but  of  a  fundamental  limitation 
due  to  the  uncertainty  principle.1 

Advances  in  theoretical  calculations  have  opened  new  doors 
in  this  area.  Calculations  do  not  have  the  same  limitations  as 
experimental  techniques.  The  properties  of  species  corresponding  to 
any  point  on  a  potential  energy  (PE)  surface  can  be  estimated,  not 
just  those  of  stable  species.  Thus,  calculations  can  estimate  the 
properties  of  transition  states.  The  word  "estimate"  is  required  in 
this  context,  however,  because  theoretical  calculations  do  have 
other  limitations . 

If  the  Schrodinger  equation  could  be  solved  exactly  for  any 
system  of  interest,  chemists  would,  in  principle,  be  able  to 
unequivocally  explain  any  chemical  process.2  For  anything  other  than 
the  hydrogen  atom  or  hydrogen  molecular  cation,  this  is  not 
possible.  As  a  result,  calculations  are  simply  one  more  tool  to  be 
used  in  concert  with  others  to  glean  information  concerning  exactly 
what  transpires  during  a  reaction.  Furthermore,  to  obtain  this 
tool,  approximations  have  to  be  made.  Different  quantum  mechanical 


models  result  from  acceptance  of  varying  sets  of  approximations.1. 
The  acceptability  of  a  model,  i.e.  a  specific  set  of 
approximations,  for  a  given  study  must  be  judged  on  the  basis  of  two 
criteria:  (1)  the  feasibility  of  applying  the  model  in  the  study, 

and  (2)  the  accuracy  of  the  model  in  predicting  properties  for  the 
types  of  species  involved  in  the  study. 

The  feasibility  of  using  a  model  is  most  commonly  determined 
by  the  time  required  for  the  calculations,  a  model  being  without 
value  if  the  calculations  take  an  unreasonable  length  of  time. 
While  some  models,  e.g.  the  Huckel  model,  can  be  applied  to  some 
problems  without  the  use  of  computers,  most  models  require  computers 
for  even  the  simplest  studies.  Even  with  computers,  the  use  of  some 
models  is  unfeasible  or  feasible  only  if  further  simplifications, 
such  as  symmetry  or  fixed  geometry,  are  imposed.  There  must  be, 
however,  substantial  evidence  from  experimental  studies  for  these 
simplifications,  or  their  effects  must  be  thoroughly  tested  for  and 
found  to  be  negligible  or  predictable  with  respect  to  the  aspects  of 
the  chemical  process  being  investigated.3  Otherwise,  they  will  be 
self-defeating,  enabling  use  of  the  model  but  casting  doubt  on  the 
conclusions  drawn  from  the  calculations. 

Even  if  its  use  in  a  particular  study  is  feasible,  a  model 
is  still  of  little  value  if  it  gives  inaccurate  results  for  the  type 
of  chemical  process  of  interest.  Since  all  models  are  based  on  sets 
of  approximations,  the  accuracy  of  results  cannot  be  assumed,  but 


must  be  established  by  comparison  of  calculated  values  of 
observables  with  experimental  values.1  The  more  properties  of  stable 
species  estimated  well  by  a  specific  model,  the  more  likely  it  is 
that  the  model  will  provide  good  descriptions  of  species 
corresponding  to  other  points  on  the  PE  surface. 

The  application  of  a  quantum  mechanical  molecular  model, 
then,  involves  three  phases:  (1)  development  of  the  set  of 
approximations  which  define  the  model,  (2)  testing  of  the  model,  and 
(3)  use  of  the  model  in  the  investigations  of  specific  chemical 
problems.  The  following  studies  cover  all  three  of  these  aspects 
with  respect  to  the  MNDO4  and  AMI5  semiempirical  models.  The  first 
study  (Chapter  2)  involves  attempts  to  parameterize  phosphorus  for 
both  normal  and  hypervalent  compounds .  The  performance  of  AMI  in 
calculating  proton  affinities  and  deprotonation  enthalpies  is 
evaluated  in  the  second  study  (Chapter  3)  to  determine  the  validity 
and  limitations  of  using  this  model  in  studies  of  reactions 
involving  protonation  and  deprotonation.  In  the  last  two  studies, 
MNDO  and  AMI  are  used  to  study  the  condensation  reaction  in 
polyketide  biosynthesis  (Chapter  4)  and  the  mass  spectral 
fragmentation  of  cis-l-nitropropene  (Chapter  5). 
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Chapter  2 


PARAMETERIZATION  OF  MNDO  AND  AMI  FOR 
NORMAL  AND  H7FERVALENT  PHOSPHORUS 

Introduction 

MINDO/3,1  MNDO,2  and  AMI3  are  semiempirical  quantum 
mechanical  molecular  models  based  on  the  Intermediate  neglect  of 
differential  overlap  (INDO)*'5  and  neglect  of  diatomic  differential 
overlap  (NDDO)  ’  approximations  to  the  Roothaan-Hall  equations , 7 
which  themselves  are  an  approximation  of  the  Hartree-Fock  (HF) 
self-consistent  field  (SCF)  equations.*  While  the  terms  INDO  and 
NDDO  refer  basically  to  the  neglect  of  various  integrals  in  the 
solution  of  the  Roothaan-Hall  equations,  the  term  "semiempirical" 
refers  to  the  replacement  of  other  integrals  by  parametric 
functions.  The  parameter  values  are  selected  to  make  the  calculated 
values  of  a  number  of  molecular  properties  for  a  wide  range  of 
molecules  agree  as  closely  as  possible  with  experimental  values. 

The  use  of  parametric  functions  with  empirically  derived 
parameters  serves  a  three -fold  purpose.  First,  by  reproducing 
"reality" ,  as  represented  by  experimental  data  for  a  large  number  of 
molecules  of  varying  types,  an  allowance  for  electron  correlation, 
which  exists  in  reality  but  is  ignored  in  the  HF/SCF  method,9  is 


implicicly  included  in  the  calculations.  Secondly,  implicit 
compensation  is  also  made  for  the  neglected  integrals  Thirdly, 

computational  time  is  greatly  reduced.  All  of  these  results  play 
essential  roles  in  attaining  the  goal  behind  the  development  of  this 
series  of  models, 1,2,3  that  is  to  provide  a  method  of  calculating 
molecular  properties  with  sufficient  accuracy  and  reliability,  at 
low  enough  cost,  to  be  of  practical  value  in  investigating  "real 
life"  chemical  problems,  particularly  with  regard  to  areas  where 
experimental  data  are  lacking  or  impossible  to  obtain  (e.g. 

transition  state  geometries).10 

Given  that  these  models  are  based  on  numerous 

approximations,  many  of  which  are  approximations  to  approximations, 

it  is  indeed  surprising  how  successful  they  are.  Among  the 

properties  for  which  MINDO/3,  MNDO  and  AMI  have  been  shown  to 

12  3 

provide  good  estimates  are  heats  of  formation,  '  ’  molecular 

123  123  23 

geometries,  ’  ’  dipole  moments,  '  ’  ionization  potentials,  ' 

electron  affinities, 11  and  proton  affinities  and  deprotonation 

enthalpies.12  For  heats  of  formation,  results  are  comparable  with  ab 

initio  calculations  with  split  valence  basis  sets.10’13 

In  general,  the  success  of  these  methods  extends  beyond 
calculations  involving  the  organic  elements  (carbon,  hydrogen, 
oxygen  and  nitrogen) .  Other  elements  for  which  one  or  more  methods 
have  been  parameterized  include  beryllium,14  boron,15  fluorine,16 
aluminum,17  silicon,16  phosphorus,10  sulfur,20  chlorine,21 
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bromine,22  tin,23  iodine,24  mercury,25  and  lead.26 

A  limitation  has  been  evident,  however,  in  the  ability  of 
these  methods  to  handle  the  so-called  hypervalent  compounds,  i.e. 
compounds  in  which  at  least  one  atom  has  an  expanded  octet.  For 
example,  the  average  absolute  error  in  the  MNDO  heats  of  formation 
for  10  trivalent  phosphorus  compounds  (P2,  PN,  P4,  PHj,  PF^ ,  PCl^, 
P(CH3)3,  P(OCH3)3,  P2H4,  P406)  is  29.0  kcal/mol,2r  while  that  for  5 
pentavalent  compounds  (PF5,  PClj,  FjPO,  CljPO,  (CH3>3PO)  is  76.9 
kcal/mol . 27  Unlike  the  individual  errors  for  the  trivalent 
compounds,  some  of  which  are  negative  while  others  are  positive,  all 
errors  for  the  pentavalent  compounds  are  positive.  The  obvious 
conclusion  from  these  results  is  that  the  bonding  parameters  do  not 
adequately  compensate  for  the  increased  electron-electron  and 
core-core  repulsion  resulting  from  bonding  five  atoms  instead  of 
three  to  the  phosphorus.  Attempts  to  correct  this  deficiency  by 
adding  d-atomic  orbitals  (d-AO)  in  the  CNDO  approximation29  to  the 
basis  set  failed  for  sulfur30  and  phosphorus.31  When  the  value 
for  sulfur  was  allowed  to  vary  to  compensate  for  contraction  of  the 
d-AOs  with  increasing  charge  on  the  hypervalent  atom,  results  were 
much  improved.  Because  of  the  functional  dependence  of  the  /3 

d 

parameter  on  atomic  charge,  however,  the  electronic  wavefunctions 
were  nonvariationally  optimized.  As  a  consequence,  full  SCF 
calculations  were  required  for  gradient  calculations  (see  below). 
The  additional  computation  time  required  due  to  this  factor,  plus 
that  due  to  expansion  of  the  basis  set  to  include  d-AOs,  made  the 


•a 
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calculations  too  slow  to  be  of  general  use.! 


Ab  initio  calculations  have  indicated  that  although  d-AOs  do 
not  participate  in  the  bonding  of  hypervalent  phosphorus  to  the 
extent  indicated  in  sp3dn  hybridization  models,  they  do  increase  the 
calculated  stability  of  hypervalent  molecules,  and  are  required  to 


accurately  describe  the  geometries  of  the  phosphine  oxides. 


32,33,34 


Specifically,  without  d-AOs  the  P-0  bond  length  of  the  hypothetical 
H3PO  molecule  was  calculated  to  be  approximately  the  same  as  that  of 
a  P-0  single  bond,  i.e.  1.6A.34  With  the  inclusion  of  d-AOs,  the 
bond  length  decreased  to  1.47A,34  which  is  close  to  the  1.48A  of 
(CH3) 3  PO  »  The  inclusion  of  d-AOs,  and  even  higher  angular 
momentum  functions,  were  also  found  to  be  essential  to  correctly 


predict  the  order  of  stability  of  P2>  P^,  and  Pg 


34,35,36,37 


MINDO/3,  MNDO  and  AMI,  however,  have  the  flexibility  of 
parametric  functions  derived  from  experimental  data,  which  the  ab 
initio  calculations,  of  course,  do  not.  As  a  consequence,  even 
without  d-AOs,  MNDO  calculates  P^  as  more  stable  than  2P2>  and  the 
P-0  bonds  in  the  phosphine  oxides  at  approximately  the  correct 
length.23  MNDO  does  calculate  cubic  Pg  to  be  68  kcal/mol  more  stable 

than  2P. ,  however.33 

4 

Despite  this  last  discrepancy,  then,  it  appears  the 
parametric  functions  in  MNDO  are  sufficient  to  compensate  for  at 
least  some  of  the  factors  for  which  d-AOs  are  required  in  ab  initio 
calculations.  Since  addition  of  d-AOs  in  the  CNDO  approximation  did 


not  give  satisfactory  results  for  both  normal  and  hypervalent 
phosphorus  compounds,31  and  addition  in  the  NDDO  approximation  would 
increase  the  number  of  two -electron  integrals  to  be  evaluated  five 
fold,31  which  would  slow  the  calculations  to  an  unacceptable  level, 
a  different  way  of  parametrically  adding  the  flexibility  required  to 
calculate  both  normal  and  hypervalent  compounds  was  sought . 

Variable  Core -Core  Repulsion  Parameter,  a 

To  determine  where  in  the  MNDO  algorithm  added  flexibility 
might  be  the  most  productive,  two  molecular  basis  sets,  one 
consisting  of  only  trivalent  phosphorus  (Pm)  compounds  and  the 
other  of  only  pentavalent  phosphorus  (PV)  compounds,  were  used  for 
independent  parameterization  trials.  Starting  with  a  set  of  values 
obtained  by  extrapolation  of  the  parameter  values  of  other  elements, 
the  normally  optimized  seven  atomic  parameters  (U  ,  U  ,  P  ,  p  , 

s»  pp  s  p 

Z^,  Z^,  a)  were  partially  optimized  according  to  the  general 

procedure  discussed  previously3  but  using  an  improved  algorithm.39 
Experimental  data  used  for  the  optimization  included  heats  of 
formation,  dipole  moments,  ionization  potentials  and  molecular 
geometries  (Tables  1  and  2) .  The  U  and  U  parameters  were 

88  pp 

relatively  insensitive  to  the  basis  sets,  changing  only  slightly 
from  initial  values  (see  below).  All  of  the  other  parameters 


Initial  P  P 


u 

■a 

-52.00 

-52.50 

-52.00 

u 

pp 

-43.50 

-43.09 

-44.06 

K 

-8.40 

-4.12 

-9.14 

-3.40 

-6.95 

-3.40 

z 

a 

2.30 

2.45 

2.30 

z 

p 

1.55 

1.70 

1.57 

a 

2.27 

2.34 

2.52 

III  V 

differed  considerably  more  between  P  and  P  values.  The  effects 

III  V 

of  the  differences  in  the  bonding  parameters  between  P  and  P  are 

twofold.  First,  s-orbital  participation  in  the  bonding  of  Pm 

III 

compounds  was  limited  by  the  lower  f}  and  the  higher  p  for  P  as 

»  p 

Y 

compared  to  P  .  This  is  consistent  with  bond  angles  close  to  90° 
for  P111  compounds  (Table  2).  Secondly,  overlap,  and,  hence, 

y 

bonding  interactions,  were  increased  for  P  compounds  by  the  lower 

Y 

and  Z^  for  P  .  This  partially  compensated  for  increased 
repulsion  due  to  more  atoms  around  phosphorus.  The  only  purpose 
served  by  the  change  in  a  was  to  decrease  the  core -core  repulsion  in 

Y 

P  compounds.  The  major  problem  with  the  current  MNDO  parameters 

lies  in  predicting  the  heats  of  formation  of  hypervalent  compounds , 

which,  as  mentioned  previously,  are  predicted  far  too  positive.  The 
III  V 

geometries  of  P  and  P  compounds,  on  the  otherhand,  are  predicted 


fairly  well . 28  Consequently,  it  appeared  a  variable  a  parameter 


might  provide  the  extra  flexibility  required  for  hypervalent 
molecules . 

One  possibility  was  to  make  a  a  function  of  the  charge  on 
phosphorus,  similar  to  the  functional  dependence  of  /S  used  for 

d 

sulfur.30  Accordingly,  a  was  optimized  for  each  of  a  series  of 
molecules  while  holding  the  other  six  parameters  fixed,  and  a  linear 
relationship  was  looked  for  between  the  optimized  a  values  and  the 
corresponding  charges  on  the  phosphorus  atom.  However,  there  was  no 
correlation. 

Since  the  need  to  add  flexibility  was  caused  by  the 

difference  in  the  number  of  atoms  bonded  to  phosphorus  in  P111  and 
v 

P  compounds,  a  correlation  was  looked  for  between  the  optimized  a 
values  and  the  sum  of  density  matrix  terms  between  the  atomic 
orbitals  on  the  phosphorus  atom  and  all  the  orbitals  on  all  other 
atoms.  Again  the  correlation  was  poor,  atoms  not  directly  bonded  to 
phosphorus  adding  too  much  to  the  sum.  When  the  density  matrix 
terms  were  multiplied  by  the  corresponding  fock  matrix  terms  to 
dampen  out  contributions  from  atoms  distal  to  the  phosphorus , 
however,  a  reasonable  correlation  existed.  Although  the  correlation 
coefficient  was  only  0.63  for  data  on  16  molecules  (Table  3)  the 
desired  trend  was  evident,  and  the  data  did  not  take  into  account 
the  effect  of  optimizing  all  parameters  with  the  new,  more  flexible 
core-core  repulsion  term. 


Based  on  these  results,  the  core-core  repulsion  parameter 

for  phosphorus,  a  ,  was  made  a  function  of  density  and  fock  matrix 
p 

terms : 


a  -  APF  +  BPF  l  f  SV  F  | 
p  b  n  v  ^  ^ 


where  /i  was  an  atomic  orbital  on  the  phosphorus  atom,  u  was  an 
atomic  orbital  on  any  other  atom  B,  and  P  and  F  were  the  density 

J  fiV  [IV 

and  fock  matrix  terms  between  n  and  v.  APF  and  BPF  were  new 
parameters  to  be  optimized. 


An  undesirable  side  effect  of  this  change  soon  became 
obvious.  For  geometry  optimization,  derivatives  of  the  energy  with 
respect  to  geometrical  parameters,  (q^),  must  be  calculated: 


aETOT  vr  3ETOT  3Vnv 

WT~ +  S  aqT" 

i  nv  pv  i 


where  E  is  the  sum  of  the  electronic  energy,  E  ,  and  the 
tot  «l 

core-core  repulsion  energy,  E  ,  and  P 1  has  it  previous  meaning. 

Equation  (2)  can  be  broken  down  into  its  components: 


3E  . 

+  S  ST- 

fiV  /ill 


dz 

+  H  er- 

fJV  /il 


/  V  v’  * 


For  variationally  optimized  electronic  wavefunctions ,  the  energy  is 
first  order  invariant  to  small  changes  in  the  density  matrix  (since 
the  energy  is  minimized  with  respect  to  the  wave  function  in  an  SCF 
cycle) ,  and  all  terms  in  the  double  summation  on  the  right  hand  side 
(RHS)  of  equation  (3)  are  zero.  E  normally  has  no  functional 

C 

dependence  on  the  density  matrix,  so  the  terms  in  the  double 
summation  on  the  RHS  of  equation  (4)  are  zero,  also.  Derivatives 
can,  therefore,  be  calculated  with  a  fixed  density  matrix. 

The  functional  variation  of  a  changed  this,  however.  While 

p 

the  electronic  wavefunction  remained  variationally  optimized,  the 
dependence  of  a  on  the  density  matrix  prevented  the  terms  in  the 
double  summation  in  equation  (4)  from  being  necessarily  zero. 
Hence,  SCF  calculations  had  to  be  used  to  calculate  the  derivatives 
of  E  .  This  made  geometry  calculations  much  slower.  This  was 
already  the  case  for  some  calculations,  e.g.  half -electron40 
calculations  of  radicals,  though,  and  the  problem  was  not  compounded 
by  an  expanded  basis  set,  as  it  was  with  sulfur.30 

Several  sets  of  parameter  values  were  used  as  initial  points 
for  optimization  of  the,  now,  eight  atomic  parameters  (U  ,  U  ,  0  , 

ss  pp  s 

0  ,  Z  ,  Z  ,  APF,  BPF)  using  the  functional  dependence  of  a  shown  in 

p  *  p  p 

equation  (1) .  The  molecules  used  in  the  basis  set  were  also  varied. 
As  previously  indicated,31  the  region  of  the  parameter  hypersurface 
in  which  the  parametrization  became  fixed  was  strongly  dependent  on 
these  selections,  as  well  as  on  the  weights  assigned  to  the  various 
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observables  used  for  the  parameterization.  While,  for  the  first 
in  v 

time,  some  P  and  P  compounds  were  calculated  relatively 
accurately  with  the  same  set  of  parameters,  totally  satisfactory 
results  were  never  obtained. 

A  major  problem  area  involved  compounds  with  multiple  bonds 
to  phosphorus,  such  as  P ^  and  PN.  Optimized  a  values  for  these 
compounds  were  surprisingly  high,  similar  to  values  required  for  PV 
compounds,  while  the  sum  of  the  product  of  density  and  fock  matrix 
terms  were  far  too  small  to  give  high  values  for  a  .  These 

p 

molecules  are  not  that  significant,  and  it  would  have  been  no  great 
loss  if  the  algorithm  for  worked  for  all  molecules  except  them. 
But  these  results  did  raise  an  interesting  question:  In  what  way 
are  and  PN  more  similar  to  PV  compound  than  Pm  compounds?  The 
answer  lies  in  the  amount  of  s -orbital  participation  in  bonding. 
Assuming  a  bonding  picture  for  P,,  and  PN  similar  to  that  for  N2 , 
i.e.  a  a  bond  of  sp  hybrids  and  two  P  -P  bonds,  s- orbital 

n  n 

participation  in  bonding  is  significant. 41  Similarly,  s-orbital 
participation  in  bonding  of  Pv  compounds  is  also  significant,  while 
that  in  normal  pm  compounds  is  slight. 

Given  these  considerations,  a  linear  correlation  was  looked 
for  between  the  optimized  a  values  used  previously  and  the  ratio  of 
interactions  with  the  phosphorus  s-orbital  to  those  with  the 
phosphorus  p-orbitals,  using  the  same  type  of  summation  as  before 
for  measuring  interactions: 
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New  Parameterization  Algorithm 

At  the  same  time  the  fatal  deficiency  in  the  calculations 
involving  the  variable  a  parameter  for  phosphorus  was  becoming 
evident,  a  new  parameterization  algorithm43  which  was  expected  to  be 
more  effective  than  the  previously  used  method,  became  available. 
The  older  method,39  which  was  used  for  all  recent  parameterization, 
did  not  involve  the  calculation  of  full  second  derivatives  for 
determining  the  search  direction  during  optimization.  As  a  result, 
the  search  directions  determined  in  the  program  were  very 
tentative.44  Parameterization  was  typically  started  by  allowing  only 
one  parameter  to  vary  at  a  time.  As  the  optimization  began  to 
stabilize,  more  than  one  parameter  would  be  allowed  to  change 
simultaneously.  Movement  around  the  parameter  surface  was 
relatively  slow,  and  could  be  influenced  by  the  researcher.  The 
decision  as  to  when  to  call  the  parameterization  complete  was 
subjective,  dependent  upon  attaining  results  which  gave  optimal 
agreement  with  experimental  data  and  with  chemical  intuition. 

In  contrast,  the  new  parameterization  included  the 
calculation  of  full  second  derivatives,  and  could,  therefore,  "see" 
the  parameter  surface  more  fully.  As  a  result,  a  more  powerful 
optimization  routine  could  be  employed  than  in  the  previous  program, 
and  the  program  was  designed  to  work  more  autonomously.  All 
parameters  set  to  optimize  were  varied  simultaneously,  and  the 
program  foui.d  the  best  mathematical  solution  (the  global  minimum) 


based  on  the  experimental  data  used  as  input.44 

When  the  new  parameterization  method  was  used  to  optimize 
all  12  MNDO  parameters  for  P'1'11  and  Pv  combined,  including  the 
one-center,  two-electron  integrals,  results  were,  at  first, 
encouraging.  The  program  clearly  reached  areas  of  the  parameter 
surface  which  the  previously  used  method  would  never  have  reached 
had  it  not  started  in  the  general  vicinity.4®  The  average  absolute 
error  in  the  heats  of  formation  for  molecules  in  the  basis  set 
typically  would  decrease  to  below  15  kcal/mol  within  at  least  20 
cycles  of  the  parameterization  (a  one -night  run) ,  even  when  the 
average  error  for  the  initial  set  of  parameters  was  over  80 
kcal/mol.  Improvements  in  ionization  potentials,  dipoles  and 
geometries  were  similar.  If  the  same  molecular  basis  set  was  used  a 
second  time  but  with  a  significantly  different  initial  set  of 
parameter  values,  the  program  would  converge  to  essentially  the  same 
set  of  values  as  previously,  indicating  It  was  the  global  minimum. 

Results  were,  however,  entirely  basis  set  dependent. 
Agreement  with  any  experimental  data  not  used  in  obtaining  a  given 
set  of  parameter  values  was  often  extremely  poor.  For  example, 
while  the  average  error  in  the  heat  of  formation  of  molecules  in  the 
basis  set  might  be  under  12  kcal/mol,  the  inversion  barrier  for  PH^ 
could  be  overestimated  by  90  kcal/mol  and  the  proton  affinity 
underestimated  by  50  kcal/mol.  If  PH+  was  then  added  to  the  basis 
set,  the  set  of  parameter  values  obtained  might  again  give  good 


results  for  the  basis  set  molecules,  including  PH^,  but  the  proton 
affinity  for  PCCH^)^  could  be  underestimated  by  40  kcal/mol. 

The  strength  of  the  new  parameterization  was,  at  once,  also 
a  weakness.  It  was  able  to  move  around  the  parameter  surface  so 
freely,  it  appeared  to  curve  fit  too  specifically  to  the  basis  set. 
Despite  many  different  approaches  to  using  the  program,  such  as 
restricting  different  parameters  from  optimizing  during  alternate 
runs  (to  try  to  stabilize  in  a  specific  region  of  the  surface)  or 
optimizing  only  P  or  P  ,  it  was  not  possible  to  influence  the 
parameterization;  once  a  basis  set  was  selected  and  weights 
assigned  to  the  observables,  the  program  was  in  complete  control. 
In  the  end,  no  satisfactory  set  of  parameter  values  was  obtained. 

Conclusions 

The  primary  shortcoming  of  MNDO  calculations  for  hypervalent 
compounds  is  the  inability  of  calculated  bonding  interactions  to 
compensate  for  the  increased  repulsion  due  to  a  larger  number  of 
atoms  around  the  hypervalent  atom  than  in  normal  valency  compounds. 
While  MNDO  has  been  shown  to  handle  pentavalent  adducts  of 
silicon,46  effectively  reproducing  the  three-center  four-electron 
bonds  thought  to  best  describe  the  primary  bonding  of  hypervalent 

32  33  34 

species,  ’  ’  there  is  no  assurance  the  stabilization  derived 
from  forming  the  adduct  is  fully  estimated  by  the  MNDO  calculations. 
The  problem  is  even  more  difficult  for  phosphorus  than  for  silicon. 
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While  all  silicon  valence  atomic  orbitals  are  involved  in  bonding  in 
normal  valency  compounds  and  in  the  pentavalent  adducts ,  this  is  not 
true  for  phosphorus.  Consequently,  parameter  values  must  ideally 
restrict  s-orbital  participation  in  the  bonding  of  pyramidal  Pm 

Y 

compounds  while  allowing  significant  bonding  participation  in  P 
compounds  and  compounds  involving  multiple  bonds  to  phosphorus . 
Unfortunately,  this  flexibility  is  difficult  to  produce  in  the  MNDO 
and  AMI  formalisms. 

To  compensate  for  these  factors,  the  core-core  repulsion 
parameter  was  modified  to  decrease  the  core-core  repulsion  energy  as 
bonding  interactions  with  phosphorus  increase.  The  difficulty, 
however,  was  finding  a  suitable  measure  of  bonding  interactions. 

The  two  functions  of  density  and  fock  matrix  terms  tried  were 
ultimately  found  to  have  fatal  flaws.  It  is  difficult  to  see  any 
other  functional  dependence  that  would  have  a  better  chance  of 
success . 


The  best  solution  still  appears  to  be  the  addition  of  d-AOs. 
While  previous  attempts  to  add  d-AOs  for  phosphorus31  and  sulfur30 
failed  to  give  satisfactory  results,  the  CNDO  approximation  under 
which  they  were  employed  was  no  doubt  a  major  part  of  the  problem. 
Ab  initio  results  indicate  exchange  terms  are  very  important  in  the 
involvement  of  d-AOs.33  Consequently,  a  better  approximation  will 
probably  be  necessary  to  incorporate  d-orbltals  effectively. 
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While  efforts  to  parameterize  for  phosphorus  using  a  new, 
"more  powerful"  parameterization  algorithm  failed  to  yield  a  useful 
set  of  parameter  values,  they  did  serve  to  point  out  the 
essentiality  of  a  chemist's  input  into  the  parameterization.  As 
pointed  out  earlier,31  rather  than  being  a  smooth  surface  with  one 
"true"  minimum,  the  parameter  surface  no  doubt  has  many  minima. 
While  they  may  all  be  stationary  points  (i.e.  mathematical 
solutions) ,  they  will  not  all  necessarily  correspond  to  solutions 
which  will  make  chemical  sense.  Consequently,  the  researcher  must 
have  the  ability  to  influence  the  parameterization.  When  a  program 
becomes  so  powerful  that  this  ability  is  restricted,  there  is  a  good 
chance  that  the  solution  obtained,  though  mathematically  correct, 
may  not  be  chemically  useful. 


TabU  1.  (Contd.) 


Molecule 

Beat  of 

a 

Formation 

Dipole 

b 

Momenta 

Ionization 

c 

Potentials 

d 

e 

p 

F  PO 

-289.5 

1.87 

13.5 

3 

j 

CH  Cl  PO 

-132.9 

3  2 

d 

i 

Cl  PO 

-133.5 

2.54 

3 

q 

f 

Br  PO 

-94.0 

11.0 

3 

1 

FF  3 

1.32 

4 

d 

f 

PF 

-381.4 

15.5 

5 

d 

f 

PCI 

-89.6 

10.7 

5 

0 

P  0 

-672 

4  10 

kcal/mol.  debye.  Vertical  Ionization  potential  in  electron 

volte  (eV).  Ref.  *7.  Ref.  *8.  Ref.  49.  ®Ref.  SO. 
h  i  j  k  1 

Ref.  SI.  Ref.  52.  Ref.  S3.  Phoephirane.  Ref.  S4. 
m  nop 

Cyclopropyl  phoaphine.  Phoaphabenzene.  Ref.  SS.  Ref.  56. 
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Table  2.  Experimental  geometrical  parameters  of  phosphorus 

a 

compounds . 


Gaomatrlcal  Parameters 


P<CH3,3 


PP  1.894 

PP  2.21 
PB  1.433 

PH  1.418;  HFH  91.7 
PH  1.420;  HPH  93.3 


1  2  12 
PP  2.219;  PH  1.414;  PH  1.417;H  PH  92.0; 

H  PP  94.3;  H  PP  99.1;  4*  74.0 


CP  1.863;  PH  1.414;  CH  1.093;  CPH  97.5; 

HPH  93.4;  HCH  109.7;  mathyl  group  tiltad 
away  from  PH^  moiaty  2.00 

CP  1.843;  CH  1.112;  CH  1.090;  PCH  111.4; 

PCH  109.8;  H  CH  108.2*  H  CH  109.4;  CPC  98.9 
a  s  a  a  a 

(a  ■  in  plana,  a  ■  out  of  plana) 


C  H  PH° 

2  4 

CP  1.867;  PH  1.428;  CC  1.502;  CH  1.092; 

cis 

CH  1.093;  CPH  95.2;  CCH  118.0; 

trans  els 

CCH  117.5;  HCH  114.4 

58 

CH  PH 

2 

CP  1.673;  PH  1.420;  CH  1.09;  CH  1.09; 

els  trans 

H  CP  124.4;  B  CP  118.4;  CPH  97.4 

els  trans 

59 

d 

C  H  P 

5  5 

2  3  3  4  2  6 

CP  1  733;  C  C  1.413;  C  C  1.384:  C  PC  101.1; 

23  234  3  4  S 

PC^C  124.4;  C  C  C  123.7;  CCC  122.8; 

PC  H  118.2;  CCH  114.7 

60 

CP 

CP  1.56 

61 

PCH 

CP  1.542;  CH  1.068;  PCH  180 

58 

PH 

PH  1.491 

58 

PO 

PO  1.474 

58 

PO 

2 

PO  1.467;  OPO  135.3 

62 

P  O 

4  S 

PO  1.638;  POP  126,4;  OPO  99.8 

60 
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Molacula 

Geometrical  Parameters 

Reference 

PHF 

2 

PH  1.412;  PF  1.582;  HPF  96.3;  FPF  99.0 

58 

PF 

3 

PF  1.563;  FPF  97.7 

58 

PF  CN 

CP  1.811;  PF  1.567;  CN  1.158;  FPF  99.1; 

2 

FPC  97.1;  PCN  171.5 

58 

PF  OCH 

PF  1.591;  PO  1.560;  CO  1.446;  CH  1.090; 

2  3 

FPF  94.8;  OPF  102.2;  COP  123.7;  HCH  110.5 

58 

pf  ra 

2  2 

PF  1.587;  PN  1.650;  NH  1.002;  NH  0. 

cia  crane 

FPF  94.6;  FPN  100.6;  PNH  123.1;  PNH 

cia  trana 

HNH  117.2 

981; 

119.7; 

58 

PF  Cl 

2 

PF  1.571;  CC1  2.030;  FPF  97.3;  FPC1  99.2 

58 

PCI  2.043;  C1PC1  100.1 

58 

PBr 

3 

PBr  2.220;  SrPBr  101.0 

60 

PI3 

PI  2.46;  IPX  100 

63 

PS 

PS  1.92 

61 

HF  PO 

2 

PH  1.387;  PF  1.539;  PO  1.437;  HPO  117.9; 

FPO  116.3;  FPF  99.8 

58 

F  PO 

3 

PO  1.437;  PF  1.522;  FPF  101.14 

58 

Cl3PO 

PO  1.455;  PCI  1.989;  C1PC1  103.7 

56 

PF 

3 

PF  1.576;  PF  1.529 

ax  aq 

64 

PCI  F 

2  3 

PF  1.539;  PF  1.591;  PCI  2.001; 

F  Sf  89.3;*C1PC1  122.2 

ax  aq 

65 

PCI 

3 

PCI  2.124;  PCI  2.020 

ax  aq 

60 

P  0 
a  10 

12  2 

PO  1.429;  PO  1.604;  PO  P  123.5; 

2  2  2  1 

0  PO  101.6;  O  PO  116.5 

60 

HF  PS 

*» 

PH  1.392;  PF  1.551;  PS  1.867;  SPF  117.4;  SPH  119.2; 

2 

SPH  119.2;  FPF  96.6 

58 
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Table  2. 

(Contd. ) 

Molecule 

Geometrical 

Parameters 

Reference 

F  PS 

PS  1.B7;  PP  1.33; 

FPF  100.3 

58 

3 

Bond  lengths  are  In  angstroms.  Bond  angles  and  dihedral  angles  are 
b 

in  degrees.  ^  Dihedral  measured  from  cis  position.  Hydrogens 

designated  H  are  the  nearly  eclipsed  pair.  Phosphirane. 
d 

Phosphabenzene . 


Table  3.  Optimized  a  values  and  corresponding  values  for  functions 

a 

of  the  density  and  fock  matrix  terms. 
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Chapter  3 


EVALUATION  OF  AMI  CALCULATED  PROTON  AFFINITIES 
AND  DEPROTONATION  ENTHALPIES 

Introduction 

Proton  transfer  reactions  play  a  basic  role  in  chemistry,  in 
particular  in  biochemistry.  A  knowledge  of  the  proton  affinities 
(PA)  of  bases,  and  of  the  deprotonation  enthalpies  (DPE)  of  acids, 
is,  therefore,  essential.  While  major  progress  has  been  made  in 
recent  years  in  the  development  of  experimental  techniques  for 
measuring  PAs,  these  have  necessarily  been  limited  to  the  gas  phase, 
and  results  are  available  only  for  a  very  limited  number  of  ions  and 
molecules1.  If  PAs  could  be  calculated  theoretically  with 
sufficient  accuracy  by  some  quantum  chemical  procedure ,  this  would 
be  of  major  value,  because  calculations,  if  feasible,  can  be  carried 
out  much  more  quickly  and  at  much  less  cost  than  experiments,  and 
they  are,  of  course,  not  limited  by  the  physical  properties  of  the 
species  being  studied. 

In  connections  such  as  this ,  the  only  theoretical  procedures 
that  need  to  be  considered  are  ab  initio  ones  based  on  the 
Roothaan-Hall2  SCF  MO  approximation  and  the  semiempirical  SCF  MO 
methods  (MINDO/33,  MNDO4 ,  AM1S)  developed  here.  Other  alternatives 


(e.g.  CNDO/2)  are  too  inaccurate  and  unreliable*. 

Numerous  ab  inicio  calculations  of  absolute  and  relative  PAs 
and  DPEs  have  been  reported  for  compounds  derived  from  carbon, 

7-9 

hydrogen,  nitrogen  and  oxygen  Many  of  these  can,  however,  be 
eliminated  from  consideration  because  of  failure  to  optimize 
geometries.  As  would  be  expected,  the  results  from  the  others 
depend  greatly  on  the  basis  set  used  and  on  the  allowance  (if  any) 
made  for  electron  correlation.  While  the  number  of  calculations  is 
too  large  for  the  results  to  be  compared  here  in  detail,  it  seems 
clear  that  a  reasonably  large  basis  set  must  be  used  if  the  average 
errors  in  DPEs  and  PAs  are  to  be  kept  within  reasonable  limits,  e.g. 
10  kcal/mol.  In  the  case  of  anionic  bases,  it  is  also  necessary  to 
include  diffuse  AOs  in  the  basis  set7 .  Calculations  at  this  level 
become  very  expensive  for  larger  molecules,  and  the  results  are 
still  by  no  means  exact. 

MNDO4  is  now  a  well  established  procedure  for  calculating 
molecular  properties10.  While  it  was  parametrized  to  reproduce 
ground  state  properties  of  neutral  closed- shell  molecules,  it  also 
gives  good  results  for  radicals,  carbenes,  and  positive  and  negative 
ions .  The  only  exceptions  are  anions  in  which  the  charge  is 
concentrated  on  a  single  atom,  where  the  calculated  heats  of 
formation  are  much  too  positive11.  This  error  is  probably11  due  to 
failure  to  allow  for  orbital  expansion  accompanying  large  localized 
negative  charges.  It  was,  therefore,  reasonable  to  hope  that  MNDO 
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might  provide  satisfactory  estimates  of  DPEs  and  PAs,  and  two  recent 
studies7,11  have  indeed  shown  that  the  results  are  comparable  with 
those  from  ab  initio  procedures  that  require  a  thousand  times  more 
computing  time.  The  errors  were,  however,  larger  than  desirable, 
particularly  in  the  case  of  certain  molecules  involving  features 
known  to  present  problems  in  MNDO.  These  problems  have  now  been 
mostly  overcome  in  a  new  "third  generation"  semiempirical  model, 
AMI5.  In  this  study,  then,  the  performance  of  AMI  in  calculations 
of  DPEs  and  PAs  is  tested  using  calculation  results  for  a  large 
number  of  molecules  for  which  apparently  satisfactory  experimental 
values  are  known. 


Procedure 

The  calculations  were  carried  out  using  the  standard  AMI5 
procedure,  as  implemented  in  the  AMPAC  package  of  computer 
programs12.  All  geometries  were  optimized  by  minimizing  the  energy 
with  respect  to  all  geometrical  variables  using  the  DFP  method13 
incorporated  in  AMPAC  and  without  making  any  assumptions .  Various 
geometries  of  protonated  methane  and  ethane  were  characterized  as 
either  true  minima,  transition  states,  or  "hilltops"  by 
diagonalizing  their  Hessian  (force  constant)  matrices  and  looking 
for  zero,  one,  or  two  or  more,  negative  eigenvalues,  respectively14. 
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DPEs  and  PAs  were  found  by  difference  from  the  calculated 
heats  of  formation  of  the  parent  molecule  and  the  derived  cation  or 
anion.  Since  AMI  gives  a  very  poor  estimate  of  the  heat  of 
formation  of  H+  (calcd.  314.9;  obsd.  367. 21 6  kcal/mol) ,  the 
experimental  value  was  used  in  calculating  DPEs  and  PAs. 

Results  and  Discussion 
A.  Proton  Affinities 

The  PA  of  a  compound  (B)  is  defined  as  minus  the  heat  of 
reaction  for  its  combination  with  a  proton; 

B  +  H+  -»  HB+ 

PA(B)  -  AHf(H+)  +  AHf(B)  -  AHf(HB+)  (1) 

Table  1  lists  calculated  and  experimental  PAs  for  60  compounds 
together  with  errors  (SAH^HB*))  in  the  heats  of  formation 
calculated  for  the  ionic  species,  these  being  found  by  using  Eq  (2); 

6AHf(HB+)  -  «AHf(B)  -  «PA(B)  (2) 

where  5AH  (B)  and  5PA(B)  are,  respectively,  the  errors  in  the  heat 
of  formation  and  in  the  PA  of  B.  Included  are  9  carbon  bases  (the 
highest  calculated  PA  is  for  protonation  of  a  carbon  atom) ,  33 
nitrogen  bases,  12  oxygen  bases,  and  6  bifunctional  bases  discussed 
separately  because  of  hydrogen  bonding  considerations. 
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Table  1.  Comparison  with  Experiment  of  AMI  Hosts  of  Formation  and  Proton  Affinities  (keal/mol). 
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Proton  Affinity 


b 

Calcd  Ezpti  Error 


Exptl 

AH  (B) 
f 


Error  in 

Calcd  AH 
_ f 

B  HB* 


Carbon  Bases 


CH 

4 


CH  CH 
3  3 


CH  »CH 
2  2 

CH  CH-CH 
3  2 


HC«CH 

CH  C*CH 
3 


•f 

CH 

3 


CH  CH* 
3  4 


CH  CH* 

3  2 

+ 

CH  CHCH 
3  3 

CH  CH  CH+ 
3  2  2 


H  C-CH+ 

2 

CH  C-CH 
3  2 


CH  CH-CH 


-8.8 

224.4 

134.0 

132.0 

2.0 

-17.8 

9.0 

7.0 

c 

134.7 

-0.7 

9.7 

a 
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a 

54.8 

261.5 

160.5 

153.3 

7.2 

54.5 

0.3 

-6.9 

a 

43.4 

233.7 

176.9 

182 

-5 

44.2 

-0.8 

4 

251.2  159.4 


i 


■WU'WIW 


WWW.’ 


um uf mi m wji »y «i ttt« ww w wwwotwtotoct1  ^ tx’h.-'v ^ c <- v* a.- ^ *.- a- ^  'it«.;.%;.-^ 


Page  44 


Tabla  1.  (Contd . ) 


Error  in 


Calcd 

AH 

£ 

Proton  Affinity 

Exptl 

Calcd 

AH 

f 

Jj 

+ 

B 

BB 

B 

HB 

Calcd 

Exptl 

Error 

AH  (B) 
f 

B 

HB 

NH 

£ 

178.9 

209.0 

(CH  )  NH 

3  2 

(CH  )  NH 

3  2  2 

a 

-5.6 

130.7 

210.9 

220.6 

-9.7 

-4.4 

-1.2 

8.5 

CH  CH  N(CH  )H 

3  2  3 

+ 

CH  CH  N(CH  )H 

3  2  3  2 

-12.6 

140.3 

214.3 

222.8 

-8.5 

f 

-11 

-2 

7 

(CH  CH  )  NH 

3  2  2 

+ 

(CH  CH  )  NH 

3  2  2  2 

-17.8 

131.5 

217.9 

225.9 

-8.0 

-17.3 

-0.5 

7.5 

H 

l 

6 

H 

+12 

N 

o 

-10.4 

139.9 

216.9 

225.2 

-8.3 

-0.8 

-9.6 

-1.3 

H 

1 

N 

O 

H 

+  12 

N 

a 

-19.0 

128.4 

219.8 

226.4 

-6.6 

-11.3 

-7.7 

-1.1 

H 

» 

(CH3)2  N  (CH3)2 

u 

H 

+*2 

(CH3 )2  N  (CH3 )2 

U 

-30.0 

109.4 

227.8 

231.7 

-3.9 

h 

-38 

8 

12 

(CH  )  N 

3  3 

+ 

(CH  )  NH 

3  3 

a 

-1.7 

152.0 

213.5 

225.1 

-11.6 

-5.7 

4.0 

15.6 

(CH  >  NC  H 

3  2  2  3 

(CH  )  N(C  H  )H 

3  2  2  3 

-6.7 

143.4 

217.1 

227.5 

-10.4 

a 

-11 

4 

14 

CH  H(C  H  ) 

3  2  3  2 

CH  N(C  H  )  H 

3  2  5  2 

-11.6 

135.0 

220.6 

230.0 

-9.4 

e 

-17 

5 

14 

CH 
l  3 

N 

o 

H  CH 
'S'  3 

6 

-6.4 

142.9 

217.9 

228.7 

-10.8 

e 

-0.5 

-5.9 

4.9 

CH 

1  3 

N 

a 

H  CH 
\  +  /  3 

0 

-15.2 

131.3 

220.5 

229.7 

-9.2 

f 

-12 

-3 

6 

»wa 

Ay*>* 

^dj 


.-.vV* 


VV, 


Page 


Tab la  1.  (Contd. ) 


Error  in 

Calcd  AH 

0 

Proton  Affinity 

Calcd  AH 

Calcd  Exptl  Error  AH  (B) 


13  13 


R«.  5.  Uniats  otherwise  notad,  experimental  FAt  tra  from  Rtf.  36«.  Rtf.  37.  Whils 
raportad  In  Rtf.  38b,  thata  vslutt  art  not  incorporsttd  Into  tba  avaluatad  seals  of  FAs  of 
Rtf.  ^36.  Uniats  otherwise  notad,  experimental  hsati  of  formation  art  from  Ref.  38.  *Ref. 
13.  Estimated  in  Ref.  36a.  ®Rtf.  39.  Ret.  36a  and  references  therein. 


(1)  Carbon  Bases.  The  average  unsigned  error  in  the  PA s 
calculated  for  the  carbon  bases  is  2.9  keal/mol16,  while  that  in  the 
heats  of  formation  of  the  conjugate  acids  is  3.9  keal/mol16.  These 
are  similar  to  the  corresponding  error  (3.1  keal/mol)  in  the  heats 
of  formation  calculated  for  the  parent  hydrocarbons . 


SCI) 

224.4(0.0) 


S(II) 

224.4(0.1) 


4v 

226.3(3.7) 


2v 

226.6(1.1) 


3b 

226.6(11.7) 


3v 

231.0(29.1) 


Figure  1.  Six  structural  wars  optimlzad  for  CB  using  AMI.  AMI  haats  of  formation,  in 
kcal/mol. ,  ara  givan  below  aach  structura.  Valuas  In  parantbases  ara  ab  lnltlo  relative 
anargias  from  Raf.  8p,  also  in  kcal/mol. 


Heats  of  formation  were  calculated  for  six  geometries  of 
CH5+  cation  (Figure  1).  In  agreement  with  high  level  ab  initio 

g_  8nxixi  17  18 

calculations  p'  '  ,  the  Cg^  structure  is  predicted  to  be 

most  stable.  AMI  calculates  the  C  , _ _ .  structure  to  be  essentially 

s(II) 

equivalent  in  stability,  its  heat  of  formation  differing  from  that 
for  by  only  0.001  kcal/mol.  While  the  AMI  results  agree  with 

the  ab  initio  calculations  as  to  the  structure  being  the  least 
stable  isomer  and  the  structure  being  the  second  least  stable 


+ 

Figure  2.  Two  structures  were  optimised  for  C  3.  using  AMI.  AMI  heats  of  formation,  in 
kcal/mol,  are  given  below  each  structure.  Values  in  parentheses  are  ah  initio  relative 
energies  from  Kef.  17,  also  in  kcal/ool. 

isomer,  the  magnitudes  of  the  differences  in  heats  of  formation 

between  these  isomers  and  C  are  significantly  less  for  AMI  than 

for  the  ab  initio  calculations.  Additionally,  the  C ^  structure 

calculated  by  AMI  is  a  very  slightly  distorted  trigonal  bipyramid, 

its  heat  of  formation  being  essentially  the  same  as  that  of  the  D3h 

structure.  The  calculated  force  constants  indicate  that  only  the 

and  structures  correspond  to  minima  on  the  potential 

energy  (PE)  surface,  each  of  the  others  having  at  least  one  negative 

eigenvalue.  Raghavachari  et  al  found  only  the  C  ,T.  structure  to  be 

s  1 1 7 

17  19 

a  true  minimum  ' 

Two  structures  were  calculated  for  C2H7+  (Figure  2).  In 
contrast  to  ab  initio  results8hh'17  and  conclusions  based  on  pulsed 
high  pressure  mass  spectrometry20,  AMI  predicts  the  Cg  structure 
alone  to  correspond  to  a  minimum.  The  structure  had  one  negative 
force  constant.  Similar  results  were  obtained  using  MINDO/38hh,zl . 


The  calculated  heats  of  formation  for  protonated  alkenes  in 
Table  1  are  for  classical  structures.  AMI  fails  to  predict  a 
nonclassical  structure  as  a  true  minimum  for  the  ethyl  cation* . 
Protonation  sites  for  asymmetric  alkenes  and  alkynes  are  predicted 
correctly,  as  are  the  preferred  para-  and  ortho-protonation  of 
toluene. 

(2)  Nitrogen  Bases.  The  average  unsigned  errors  in  the  PAs 
calculated  for  neutral  nitrogen  bases ,  and  for  the  heats  of 
formation  calculated  for  them  and  for  their  conjugate  acids,  are  5.8 
kcal/mol,  3.9  kcal/mol  and  5.3  kcal/mol,  respectively.  The  errors 
in  the  calculated  PAs  and  in  the  heats  of  formation  calculated  for 
the  conjugate  acids  are  again  similar  to  that  for  the  neutral  bases. 

The  conjugate  acids  for  all  amines  except  aniline  are 
calculated  to  be  unstable  relative  to  the  corresponding  bases;  i.e. 
5AH^(HB+) -SAH^(B)>0.  Furthermore,  these  differences,  and  the  errors 
in  the  PAs,  generally  increase  in  the  sequence  1°  <2°  <3°  .  The 
former  trend  is  probably  due  to  AMI  being  parameterized  for  ground 
state  neutral  molecules,  in  each  of  which  nitrogen  has  a  localized 
lone  pair  of  electrons .  The  latter  trend  can  be  attributed  to  the 
fact  that  ammonium  ions  have  four  atoms  bonded  to  nitrogen,  whereas 
the  molecules  included  in  the  parameterization  basis  set  had  a 
maximum  of  three.  Moreover,  few  of  the  latter  had  more  than  one 
alkyl  group  bonded  to  nitrogen. 


As  a  result  of  these  factors,  the  PAs  of  amines  are 
underestimated  by  AMI  in  an  irregular  manner,  much  as  they  were  by 
MNDO11.  The  errors  are  less  with  AMI  than  with  MNDO,  however,  and 
some  relative  basicities  predicted  incorrectly  by  MNDO  are  now 
predicted  correctly.  For  example,  AMI  now  correctly  predicts 
trimethylamine  to  be  more  basic  than  methylamine,  though 
dimethylamine  is  still  predicted  to  be  less  basic  than  methylamine. 

AMI  predicts  preferential  N-protonation  of  aniline,  in 
agreement  with  experiment22  and  theory  .  AMI  results  disagree  with 

8e 

the  results  of  Del  Bene's  high  level  ab  initio  calculations  , 
however,  in  predicting  cis-ethanimine  to  be  more  stable  than  the 
trans-  isomer. 

AMI  correctly  predicts  quinuclidine  to  be  more  basic  than 
its  unsaturated  analog.  The  difference  between  the  two  calculated 
PAs  is  significantly  less  than  the  experimental  results,  however, 
indicating  that  AMI  underestimates  the  inductive  effect  of  the 
sp2 -carbon  in  the  unsaturated  compound. 

(3)  Oxygen  Bases.  The  average  unsigned  error  in  the 
calculated  PAs  of  the  oxygen  bases  is  6.8  kcal/mol16,  while  that  in 
the  calculated  heats  of  formation  of  the  cations  is  5.4  kcal/mol1 s, 
the  corresponding  error  for  the  neutral  bases  being  5.1  kcal/mol. 
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Experimental  trends  within  each  group  of  compounds  (i.e. 
alcohols,  ethers  and  carbonyl  compounds)  are  accurately  reproduced 
by  AMI.  As  a  whole,  however,  alcohols  and  ethers  are  predicted  to 
be  weaker  bases  than  experiment  indicates  by  approximately  10 
kcal/mol.  As  in  the  case  of  ethanimine,  AMI  again  differs  from  Del 
Bene's  ab  initio  results8*  in  predicting  preferential 
cis-protonation  of  acetaldehyde  and  formamide. 

(A)  Bifunctional  Bases .  The  unsigned  average  error  in  PAs 
for  the  six  bifunctional  bases  listed  is  9.3  kcal/mol,  compared  with 
10.8  and  2.2  kcal/mol,  respectively,  for  the  heats  of  formation  of 
the  bases  and  conjugate  acids.  Thus,  for  these  compounds,  the  heats 
of  formation  of  the  conjugate  acids  of  the  bifunctional  bases  are 
calculated  more  accurately  than  those  of  the  neutral  bases .  Note 
that  the  errors  for  the  neutral  bases  are  all  negative  and  that  the 
errors  for  three  of  four  n-alkyl  bifunctional  bases  are 
exceptionally  large  (see  below). 

As  mentioned  previously,  these  bifunctional  bases  are  of 
particular  interest  in  that  the  PAs  for  five  of  the  six  compounds 
are  much  larger  than  those  of  alkyl  amines  of  comparable 
polarizability,  the  exception  being  hydrazine.  The  high  PAs  can  be 
attributed  to  intramolecular  hydrogen  bonding  in  the  conjugate 
acids1*’24.  For  hydrogen  bonding  to  play  a  role  in  the  protonation 
of  hydrazine,  the  conjugate  acid  would  have  to  have  a  structure 
analogous  to  the  nonclassical  structure  of  the  ethyl  cation.  No 


such  structure  was  found  to  be  a  minimus)  on  the  AMI  PE  surface. 

Two  geometries  were  calculated  for  each  of  the  four  n- alkyl 
bifunctional  bases  listed  in  Table  1,  one  a  "cyclic"  structure 
(chain  dihedral  angles  -0°)  and  the  other  an  "extended"  structure 
(chain  dihedral  angles  -180°).  It  was  assumed  that  hydrogen  bonding 
would  play  a  role  only  in  the  cyclic  structure.  Heats  of  formation 
are  given  for  both  optimized  geometries.  PAs  are  calculated  using 
the  heats  of  formation  of  the  more  stable  conformers. 

In  the  more  stable  conformer  of  1 , 2-diaminoethane ,  one  amino 
group  is  gauche  to  the  other.  However,  the  relative  orientation  of 
the  two  amino  groups  in  the  optimized  structure  precludes  hydrogen 
bonding  between  them.  The  lower  heat  of  formation  calculated  for 
this  conformer  relative  to  that  for  the  conformer  in  which  the  two 
amino  groups  are  anti  to  each  other  is  a  manifestation  of  the  gauche 
effect25’26. 

If  no  hydrogen  bonding  is  involved  in  the  protonated  base, 
the  difference  between  the  heats  of  formation  calculated  for  the 
cyclic  and  extended  conformers  of  the  acid  should  be  approximately 
the  same  as  for  the  neutral  base,  i.e.  AH^  of  the  cyclic  structure 
of  ^NCHjC^NHj  should  be  140.5  kcal/raol.  Consequently,  AMI 
calculates  hydrogen  bonding  to  stabilize  the  cyclic  structure  of  the 
conjugate  acid  by  4.3  kcal/mol. 


Table  2.  Comparison  with  Experiment  o£  AMI  Hydrogen  Bonding  Stabilization  in 
the  Conjugate  Acid*  of  n-Alkyl  B1 functional  Bases  with  Experimental  and  Ab 
Initio  Results  (keal/mol). 


e 


Conjugate  Acid 

AMI 

Exptl 

Ab  Initio 

+ 

H  N(CH  )  NH 

2  2  2  3 

4.3 

•  b  c 

6.7  ,9.1  ,12.6 

f  8 

24.8  ,7.7 

H  NCCH  )  NH+ 

2  2  3  3 

9.3 

e  b  c 

14.2  ,14.2  ,20.5 

HO(CH  )  NH+ 

2  3  3 

d 

7.7 

a  b 

8.8  ,15.2 

+ 

H  NCCH  )  NH 

2  2  4  3 

14.2 

a  b 

17.9  ,19.8 

«  b 

Raf.  28.  Ref.  29. 

°Raf .  30. 

d 

Hydrogen  bonding  stabilization  in  this 

acid  is  underestimated  due  to  hydrogen  bonding  interaction  in  the  neutral  base. 

Sm  text.  *R«r.  8jjj.  STO-  3G//STO-3G.  *4-31G/7STO-3G. 

The  stabilization  due  to  hydrogen  bonding  in  the  conjugate 

acids  of  1,3-diaminopropane,  3*aminopropan-l-ol  and 

1,4-diaminobutane  was  calculated  in  the  same  manner  as  for 

1,2-diaminoethane.  The  AMI  results  are  given  in  Table  2,  along  with 

experimental  estimates  by  Meot-Ner  et  al2*,  Buschek  et  al28,  and 

Yamdagni  et  al30 ,  and  the  ab  initio  results  of  Houriet  et  al8J  . 

The  AMI  values  are  similar  to  the  most  recent  experimental  values, 

i.e.  those  of  Meot-Ner  ’  ,  the  N-H*  •  -N  hydrogen  bond  being 

somewhat  underestimated  by  AMI.  The  calculated  value  for 

3-aminopropan-l-ol  does  not  indicate  the  full  extent  of  hydrogen 

bonding  in  the  acid  since  hydrogen  bonding  is  also  involved  in  the 
29  32 

neutral  base  ’  .  The  preferred  N-protonation  as  opposed  to 

O-protonation  is  in  agreement  with  previous  theoretical  calculations 
and  experimental  results8JJJ . 
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Tabla  3.  Comparison  of  Errors  in  AMI  Haats  of  Formation  for  n-Alkanas, 
n-Alkyl  Aminas  and  n- Alkyl  Alcohols  (kcal/mol) . 


CH  <CH  ) 
3  2 

CB 
n  3 

CH  (CH  ) 
3  2 

NH 

n  2 

CH  (CH  ) 
3  2 

OH 

n 

AMI 

Exptl 

AMI 

Exptl 

AMI 

Exptl 

AH 

f 

b 

AB 

£ 

Error 

AH 

f 

b 

AB 

£ 

Error 

AH 

f 

b 

AH 

f 

Error 

n-0 

a 

-17.4 

-20.0 

2.6 

a 

-7.4 

-5.5 

-1.9 

a 

-57.0 

-48.2 

-8.8 

n-1 

a 

-24.3 

-25.0 

0.7 

a 

-15.1 

-11.3 

-3.8 

a 

-62.7 

-56.2 

-6.5 

n-2 

a 

-31.1 

-30.0 

-i.i 

a 

-22.1 

-16.8 

-5.3 

-70. 6* 

-61.0 

-9.6 

n»3 

a 

-37.9 

-35.1 

-2.8 

-28.9 

-22.0 

-6.9 

-77.3 

-65.7 

-11.6 

n-4 

-44.1 

-39.9 

-4.2 

-35.1 

c 

-26 

-9 

-84.1 

-70.4 

-13.7 

a  b  c 

Raf.  5.  Raf.  38.  Estimatad  in  Raf.  36a. 


Despite  the  good  predictions  of  the  strengths  of  the 
hydrogen  bonds,  the  errors  in  the  calculated  PAs  for  the  n-alkyl 
bifunctional  bases  are  exceptionally  large.  The  errors  in  the 
corresponding  PAs  are,  moreover,  due  largely  to  errors  in  the 
calculated  heats  of  formation  of  the  neutral  bases  rather  than  in 
those  of  the  conjugate  acids.  The  reason  for  this  is  a  combination 
of  systematic  errors,  one  reflecting  increasing  alkyl  chain  lengths 
and  the  other  substitution  of  nitrogen  or  oxygen  for  a  primary 
carbon  atom.  Table  3  lists  the  errors  in  AMI  heats  of  formation  for 
a  series  of  n-alkanes,  n-alkyl  amines  and  n-alkyl  alcohols.  The 
errors  for  the  alkanes  change  by  an  average  of  -1.7  kcal/mol  for 
each  additional  methylene  group  in  the  alkyl  chain33.  Likewise 
substitution  of  a  nitrogen  atom  for  a  primary  carbon  changes  the 


Table  4.  Comparison  of  Ab  Initio  and  AMI 
Proton  Affinities9  (kcal/mol). 


6-31G*/ /3-21G 

AMI 

Moloculo 

a 

PA 

b 

Error 

c 

Error 

NH 

3 

218.4 

14.4 

5.3 

CH  NH 

3  2 

228.3 

14.4 

-3.0 

(CH  )  NH 

3  2 

233.2 

14.6 

-9.7 

acs 

178.1 

6.7 

12.0 

CH  CN 

3 

194.0 

3.6 

2.0 

H  0 

2 

173.3 

d 

6.9 

d 

-2.0 

CH  OB 

3 

189.3 

7.6 

-10.0 

(CH  )  0 

3  2 

198.6 

6.5 

-14.7 

H  CO 

2 

181.4 

9.7 

2.7 

CH  CHO 

3 

194.0 

7.4 

-1.7 

Avg.  Hnclgnad  Error 

9.4 

6.3 

a  b 

Ref.  8jj.  Errors  are  calculated  using 

c 

experimental  values  from  Table  1.  From 
Table  1.  Based  on  the  average  of  the  two 
experimental  values  given  in  Table  1. 

Is  probably  0.4  -  2.9  kcal/mol  higher  than  that  for  ring 

protonation.  The  small  difference  between  N-protonation  and  ring 
protonation  seems  quite  reasonable,  given  that  aniline  undergoes 
N-protonation22  in  the  gas  phase  while  m-diaminobenzeneZ2a,c  and 
1-aminonaphthalene24  protonate  in  the  ring. 


(5)  Ab  Initio  Calculations.  The  results  of  high  level  ab 
initio  calculations8,3,5  of  PAs  for  ten  of  the  molecules  listed  in 
Table  1  are  summarized  in  Table  4.  While  the  errors  naturally 
fluctuate,  the  AMI  values  are,  on  average,  more  accurate  by  a 
significant  margin. 


B.  Deprotonation  Enthalpies 

The  deprotonation  enthalpy  (DPE)  of  compound  HB  is  the  heat 
of  reaction  for  loss  of  a  proton  to  form  the  conjugate  base; 

HB  -  H+  +  B* 

DPE(HB)  -  AHf(H+)  +  AHf(B*)  -  AHf(HB)  (3) 

The  DPE  of  a  compound  is  thus  equal  to  the  PA  of  its  conjugate  base. 
Which  terminology  is  chosen  (PA  or  DPE)  depends  on  the  charges 
present,  PA  and  DPE  being  regarded  primarily  as  properties  of 
neutral  molecules. 

Calculated  DPEs  for  80  compounds  are  given  in  Table  5 ,  along 
with  corresponding  experimental  values.  Also  included  are  the 
errors  in  the  calculated  heats  of  formation  of  the  conjugate  bases, 
5AHf(B  ),  calculated  using  the  equation; 


(4) 


SAH  (B  )  -  «DPE(HB)  +  5 AH  (HB) 


where  5DPE(HB)  and  5AH£(HB)  are  the  errors  in  the  calculated  DPE  and 
heat  of  formation  of  HB,  respectively.  The  listed  compounds  include 
46  carbon  acids  (the  lowest  DPE  is  for  a  hydrogen  bonded  to  a 
carbon),  10  nitrogen  acids  and  24  oxygen  acids. 

(1)  Carbon  Acids.  The  average  unsigned  error  in  calculated 
DPEs  for  the  carbon  acids  is  8.0  kcal/mol16.  The  values  for  CH^, 
C^,  HCN,  and  the  alkynes  are  exceptionally  poor  as  a  result  of 
large  errors  in  the  calculated  heats  of  formation  of  the 
corresponding  anions.  When  these  eight  compounds  are  excluded  from 
the  statistical  analysis,  the  average  unsigned  error  is  reduced  to 
5.7  kcal/mol1®.  The  average  unsigned  error  in  the  heats  of 
formation  of  the  38  anions  is  also  5.7  kcal/mol18.  These  errors 
compare  favorably  with  the  error  of  5.1  kcal/mol  for  the  38 
corresponding  neutral  acids. 

The  poor  results  for  the  CH^  and  CN  anions  can  be 
attributed  to  the  failure  of  AMI  to  allow  for  orbital  expansion  on 
atoms  bearing  large  negative  charges8.  Similar  problems  were 
encountered  in  calculations  for  anions  using  MNDO11.  The  calculated 
heats  of  formation  are  expected  to  be,  and  are,  too  positive 
whenever  the  formal  charge  in  an  anion  is  largely  concentrated  on  a 
single  atom.  This  is  also  the  case  for  HO  5  and  H£N  anions  (see 
below).  The  same  problem  arises,  as  expected,  in  the  case  of 
allenyl  anion40  and  propynyl  anion.  Both  are  incorrectly41 
predicted  to  be  less  stable  than  propargyl  anion  (  CH0OCH,  C„  ). 


Page  61 


While  the  negative  charge  in  the  allenyl  or  propynyl  anion  is 
concentrated  largely  on  a  single  carbon  atom,  in  propargyl  it  is 
dispersed  over  an  allylic  system. 

Experiment43  indicates  that  the  propynyl  anion  is  lower  in 
energy  than  allenyl  anion  by  3  kcal/mol  while  ab  initio 
calculations8"’ 44,45  imply  the  allenyl  anion  is  lower  than  propargyl 
by  7-9  kcal/mol.  Thus  the  difference  between  the  propargyl  and 
propynyl  anions  is  ca  11  kcal/mol.  Using  the  experimental42  value 
for  the  heat  of  formation  for  the  latter,  we  arrive  at  an  estimate 
(68  kcal/mol)  for  that  of  propargyl  anion,  in  good  agreement  with 
AMI  (64.5  kcal/mol).  This  supports  the  interpretation  given  above 
of  the  errors  for  allenyl  and  propynyl. 

A  similar  situation  is  expected  in  the  case  of  propene  and 
its  derivatives.  While  allyl  anion  is  correctly  predicted  to  be  the 
most  stable  isomer,  the  difference  between  it  and  1-  or  2-propenyl 
anion  is  probably  overestimated  by  AMI.  The  relative  order  of 
stability  predicted  by  AMI  agrees  with  that  given  by  MNDO11  and  by 
the  ab  initio  calculations  of  Boerth  and  Streitwieser8™1,  except  in 
that  the  E-  and  Z-  isomers  of  1-propenyl  anion  are  predicted  by  AMI 
to  have  essentially  identical  energies .  AMI  also  reproduces  the 
effects  of  substituents  at  the  central  carbon  atom,  attributed  by 
Bartmess  and  Burnham  to  be  primarily  polar  in  nature48. 
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Table  5.  Comparison  with  Experiment  of  AMI  Boats  of  Formation  and  Dapxotonation  Enthalpies 
(kcal/mol) . 
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Raf.  5.  Unless  otherwise  noted,  experimental  DPEs  are  from  Ref.  52.  Unless  otherwise 

d  e  f 

noted,  experimental  heats  of  formation  are  from  Ref.  38.  Ref.  54.  Ref.  40.  Calculated 

8  h 

for  the  formation  of  propyne  upon  protonation.  Estimated  in  Ref.  52.  Ref.  36a  and 

references  therein.  Ref.  15  ^Ref.  lb.  See  text.  Estimated  by  group  additivit)  method, 

Ref.  S3.  Value  for  C-(H)  <C0)(0)  approximated  as  that  for  C-(H)  (0)(C).  Ref.  52  and 
m2  2 

references  therein.  Ref.  39. 


The  errors  in  the  DPEs  calculated  for  the  nitro  alkanes  are 
also  relatively  large.  These,  however,  are  due  primarily,  to  errors 
in  the  heats  of  formation  calculated  for  the  neutral  acids8  rather 
than  in  those  for  their  conjugate  bases. 


In  all  but  one  case  AMI  overestimates  the  stabilizing  effect 
of  a  methyl  or  phenyl  substituent  at  the  anionic  center  in  a 
carbanion,  the  exception  being  the  2-nitro-2-propyl  anion.  The 
effect  is  most  significant  when  the  parent  anion  is  primary,  the 
error  being  approximately  -8  keal/mol  for  methyl  substitution  and 
-10  keal/mol  for  phenyl  substitution47. 


The  experimental  evidence48  and  ab  initio7  calculations 
agree  with  the  AMI  prediction,  that  the  a-hydrogens  of  carbonyl 
compounds  are  more  acidic  than  the  acyl  hydrogen.  The  difference 
between  the  corresponding  DPEs  in  acetaldehyde  is  calculated  by  AMI 
to  be  23  kcal/mol,  in  reasonable  agreement  with  the  value  (29 
kcal/mol)  derived  from  4-31  +  G//4-31  +  G  and  MP2/4-31  +  G//4-31  +  G 
calculations7.  Additionally,  AMI  predicts  the  lowest  DPE  in 
acetonyl  acetate  to  be  that  for  the  methene  hydrogens ,  in  agreement 
with  an  earlier  assignment11*.  However,  due  to  uncertainty 
concerning  the  change  in  entropy  accompanying  the  loss  of  a  methene 
hydrogen,  neither  the  deprotonation  site  nor  the  experimental  value 
of  the  DPE  is  firmly  established49. 

(2)  Nitrogen  Acids.  The  average  unsigned  error  in  the  DPEs 
calculated  for  the  nitrogen  acids  listed  in  Table  5  is  6.1  kcal/mol. 
If  NH^  is  eliminated  for  the  reasons  indicated  above,  the  error  is 
reduced  to  4.1  kcal/mol,  and  that  for  the  heats  of  formation  of  the 
anions  to  4.6  kcal/mol.  These  values  again  compare  favorably  with 
the  corresponding  error  (4.1  kcal/mol)  in  the  heats  of  formation  of 
the  neutral  acids. 

The  relative  acidities  of  methylamine  and  ethylamine  are 
predicted  incorrectly  by  AMI.  This  error  is  due  to  two  systematic 
errors  discussed  above  in  detail,  one  being  the  error  accompanying 
methyl  substitution  at  an  anionic  site  while  the  other  involves 
increasing  alkyl  chain  length.  The  prediction  that  methylamine 


behaves  as  a  nitrogen  acid  rather  than  a  carbon  acid  coincides  with 
chemical  intuition  and  with  the  ab  initio  results  of  Lohr  and 
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The  acidities  of  p-methyl-  and  p-methoxyaniline  are 
predicted  incorrectly  relative  to  aniline.  These  errors  can  again 
be  attributed  to  the  overestimation  of  the  stabilizing  effects  of 
methyl  groups  attached  to  negatively  charged  carbon  atoms. 

(3)  Oxygen  Acids.  For  the  oxygen  acids  listed  in  Table  1, 
the  average  unsigned  error  in  the  DPEs  is  9.4  kcal/mol16.  The  error 
in  the  case  of  H^O  has  already  been  discussed.  The  errors  for 
oximes  are  also  significantly  larger  than  those  for  the  other 
compounds  (see  below).  Excluding  water  and  the  oximes  from  the 
statistical  analysis,  the  average  unsigned  error  for  the  remaining 
19  DPEs  falls  to  7.6  kcal/mol1*,  and  that  in  the  heats  of  formation 
of  the  conjugate  bases  to  7.4  kcal/mol18.  These  are  larger  than  the 
corresponding  errors  for  the  carbon  acids  or  nitrogen  acids,  and 
they  are  also  larger  than  the  average  unsigned  error  (5.0  kcal/mol) 
in  the  heats  of  formation  of  the  parent  acids.  Since  the  errors 
associated  with  compounds  containing  oxygen  are,  in  general,  larger 
than  those  for  hydrocarbons  or  nitrogen-containing  compounds5,  and 
since  the  charge  on  oxygen  in  an  oxygen-containing  anion  is  usually 
large,  it  is  not  surprising  that  the  AMI  errors  for  such  species  are 
larger  than  usual. 


I 

? 
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The  large  errors  in  the  DPEs  of  oximes  are  due  to  the 
overestimation  of  the  stability  of  the  conjugate  bases.  Given  the 
tendency  of  AMI  to  overestimate  the  stabilization  associated  with 
phenyl  and  methyl  substitution  on  the  anionic  center  (see  above) , 
the  pattern  of  errors  in  the  heats  of  formation  of  the  variously 
substituted  oxime  conjugate  bases  is  consistent  with  significant 
charge  buildup  on  the  methene  carbon  (i.e.  errors  in  heats  of 
formation  are  large  and  negative  with  methyl  and  phenyl  substitution 
on  the  methene  carbon,  but  decrease  significantly  with  t -butyl 
substitution) .  Analysis  of  atomic  charges  in  ethanal  oxime 
indicates  that  0.40  of  the  formal  charge  of  the  anion  is  on  the 
methene  carbon,  with  the  methyl  group  absorbing  another  0.13. 
Additionally,  the  coefficient  of  the  methene  carbon  p-orbital  in  the 
HOMO  increases  from  0.57  in  the  neutral  acid  to  0.75  in  the  anion. 
These  results  point  to  an  overestimation  of  charge  delocalization  by 
AMI  through  a  conjugated  system  terminated  by  a  methyl  or  phenyl 
group,  similar  to  the  effect  noted  with  the  anilines. 

While  the  decrease  in  the  DPEs  of  aliphatic  alcohols  with 
increasing  size  of  the  alkyl  group  is  reproduced  qualitatively  by 
AMI,  the  calculated  difference  (7.2  kcal/mol)  between  the  highest 
DPE  (CH30H)  and  the  lowest  DPE  (t-Bu2CH0H)  is  less  than  that 
observed  (11.9-14.1  kcal/mol).  Moreover,  the  order  of  decreasing 
DPEs  is  not  in  agreement  with  experiment.  While  the  scatter  in  the 
errors  of  the  calculated  heats  of  formation  of  the  neutral  alcohols 
may  be  partly  responsible,  the  major  problem  seems  to  lie  in  the 


errors  increasing  with  the  size  of  the  alkyl  groups.  Thus  AMI  does 
not  fully  reproduce  the  charge -induced  dipole  stabilization  of  the 
anions  to  which  the  trend  in  gas -phase  acidities  in  alcohols  has 
been  attributed30,51. 

The  benzylic  hydrogens  are  calculated  to  be  more  acidic  than 
the  hydroxyl  hydrogen  in  benzyl  alcohol.  Applying  a  10  kcal/mol 
correction  to  the  heat  of  formation  of  the  carbon  anion  to 
compensate  for  the  overestimation  by  AMI  of  the  added  stability  due 
to  phenyl  substitution  at  an  anionic  center,  the  hydroxyl  hydrogen 
becomes  the  more  acidic,  as  is  observed42. 

(A)  Comparison  with  Ab  Initio  Results.  Table  6  lists  ab 
initio  DPEs,  without  (4-31  +  G//4-31  +  G)  and  with  (MP2/A-31  + 
G//4-31  +  G)  allowance  for  electron  correlation  by  second-order 
Moller-Plesset  (MP)  perturbation  theory7,  for  nine  of  the  compounds 
included  in  Table  1.  Five  other  molecules  included  in  both  studies 
(CH^,  NHj.  H^O,  HCN,  HCCH)  have  not  been  included  in  Table  6  since 
it  has  already  been  recognized  that  AMI  cannot  describe  the 
corresponding  conjugate  bases  adequately. 

The  basis  set  in  the  ab  initio  study  includes  diffuse  AOs, 
known7  to  be  essential  for  the  proper  description  of  anions.  While 
a  few  ab  initio  values  are  better  than  the  corresponding  AMI  ones, 
the  overall  average  error  is  less  for  AMI. 


Table  S.  Comparison  of  A b  Initio  and  AMI  Daprotonation  Enthalpies9 
(kcal/mol) . 


4-31+G//4-31+G 

MP2/4-31+G//4-31+G 

AMI 

Molsculs 

DPE* 

b 

Error 

a 

DPE 

b 

Error 

c 

Error 

C  B 

2  6 

439.1 

18.1 

432.0 

11.0 

-1.9 

C  H 

2  4 

423.8 

17.8 

417.9 

11.9 

12.5 

CH  CHCH 

2  3 

403.3 

14.8 

399.1 

8.4 

-2.5 

CH  CH 

3 

386.1 

14.0 

383.4 

11.3 

6.6 

CH  CHO 

3 

374.3 

8.1 

369.2 

2.8 

5.2 

CH  NH 

3  2 

422.6 

19.4 

410.7 

7.5 

4.5 

CH  OH 

3 

393.1 

d 

12.8 

381.4 

d 

1.1 

( 

5.1 

CH  CH  OH 

3  2 

391.7 

14.3d 

378.7 

d 

1.5 

< 

6.9 

HCOOH 

346.9 

1.7 

337.3 

-7.9 

9.4 

Avg.  Unsigned  Error 

13.3 

7.0 

6.0 

a  b 

Ref.  7.  Errors  are  calculated  using  experimental  values  from  Table 
c  d 

3.  From  Table  3.  Based  on  the  average  of  the  two  experimental 
values  quoted  in  Table  3. 


Conclusions 


With  a  few  exceptions,  AMI  seems  to  be  an  effective  method 
for  studying  processes  involving  deprotonation  or  protonation  of 
neutral  molecules .  The  errors  in  the  calculated  DPEs  and  PAs ,  as 
well  as  in  the  calculated  heats  of  formation  of  deprotonated  and 
protonated  species,  are  comparable  with  those  in  the  calculated 
heats  of  formation  of  the  neutral  precursors .  Intramolecular 
hydrogen  bonding  in  protonated  bifunctional  bases  is  also 
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effectively  reproduced. 


The  main  problems  encountered  with  AMI  involve  small  anions , 
in  which  the  charge  is  largely  concentrated  on  one  atom,  and  anions 
formed  by  the  deprotonation  of  oximes .  Systematic  errors  accompany 
introduction  of  methyl  or  phenyl  substituents  at  anionic  centers . 
Systematic  errors  also  arise  in  the  extension  of  alkyl  chains  by 
addition  of  methene  groups  and  in  substitution  of  amine  and  hydroxyl 
groups  for  methyl  groups  bonded  to  secondary  or  tertiary  carbons  in 
neutral  molecules.  The  errors  involved  in  the  deprotonation  of 
alcohols  and  protonation  of  amines  are  not,  however,  totally 
systematic.  As  a  result,  the  relative  DPEs  of  alcohols  and  PAs  of 
amines  are  not  accurately  reproduced  by  AMI. 


The  comparisons  in  Tables  4  and  6  suggest  that  AMI  performs 
as  well  here  as  do  quite  high  level  ab  initio  procedures55 .  The 
fact  that  the  comparisons  refer  only  to  a  few  simple  cases  is  due 
simply  to  the  dearth  of  relevant  ab  initio  data.  Calculations  at 
this  level,  if  carried  out  with  full  geometry  optimization,  become 
very  expensive  for  larger  molecules.  Since  the  accuracy  achieved  by 
AMI  is  sufficient  for  the  results  to  be  chemically  useful  and  since 
it  can  be  used  to  study  reactions  of  quite  large  molecules  at 
moderate  cost,  AMI  should  prove  useful  as  an  aid  in  interpretating 
proton  transfer  in  chemistry  and  biochemistry. 
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Chapter  4 


An  AMI  and  MNDO  Study  of  the  Condensation  Reaction 
of  Polyketide  Biosynthesis 

Introduction 

Polyketide  biosynthesis  is  responsible  for  a  diverse  group 
of  natural  products .  Defined  in  biosynthetic  terms  as  opposed  to 
structural  terms,  polyketides,  in  general,  comprise  structures 
derived  essentially  from  poly-/?-ketomethylene  chains f CHRCO]  - ,  R 
commonly,  but  not  exclusively,  hydrogen.1  As  such,  this  category  of 
compounds  includes  fatty  acids,  as  well  as  compounds  with  more 
exotic  structures,  such  as  6-methylsalicylic  acid.2 

While  the  ultimate  products  vary  greatly,  the  biosynthesis 
of  fatty  acids  and  other  polyketides  follow  an  essentially  similar 
pathway3,4  involving  sequential  condensation  of  two-carbon  subunits 
onto  a  lengthening  carbon  skeleton.  The  diversity  of  the  products 
arises  from  variations  in  the  reactions  occurring  between  successive 
condensations . 

The  two-carbon  subunits  co.ld  be  derived  directly  from 


simple  aliphatic  acids : 
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-CH(CH3)CO  -  'propionate' 

-CH(CzH5)CO-  'butyrate' 

In  actuality,  they  enter  the  synthetic  process  as  malonic  acid 
derivatives  (HOOCCHRCOOH) ,  which  are  more  effective  nucleophiles.1 
While  many  natural  polyketides  can  be  described  as  "wholly 
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acetate - derived" ,  having  been  synthesized  from  an  initial  acetyl 
"primer"  unit  by  successive  condensations  of  malonate  groups  (which 
are  synthesized  by  carboxylation  of  acetyl  units  )  with  coincidental 
evolution  of  C02>  some  enzymes  are  known  to  accept  different  acyl 
"primer"  units  and  substituted  malonate  groups.1 

With  the  exception  of  the  transacylase  activities ,  which 
transfer  acetyl  and  malonyl  groups  between  coenzyme  A  (CoA)  and  the 
enzymes,  thiol  esters  are  used  in  place  of  alcohol  esters.5  Acyl  and 
malonyl  groups  are  stronger  nucleophiles  and  electrophiles  in  thiol 
esters  than  in  alcohol  esters .  According  to  the  resonance 
explanation  normally  given  for  this  difference  in  behavior,*’6 
nonbonding  electrons  of  the  alcohol  oxygen  atom  of  alcohol  esters 
are  delocalized  into  the  carbonyl  group,  and  all  three  resonance 
structures  shown  below  are  important.  In  thiol  esters,  the  overlap 
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of  the  3p  orbitals  with  the  it  system  of  the  carbonyl  group  is  less 
effective,  and  resonance  structure  (III)  is  less  important. 
Consequently,  the  carbonyl  carbon  is  more  electrophilic  in  thiol 
esters  than  alcohol  esters,  and  the  thiol  esters  are  better 
acylating  agents.  While  there  is  some  evidence  against  this 
argument,7  the  increased  activity  in  thiol  esters  is  not  disputed. 
Indeed,  this  increased  activity  was  evident  in  the  results  of  a  gas 
phase  study  by  Bartmess  et  ai , 8  in  which  acetaldehyde  enolate 
reacted  with  methyl  thiolformate  to  form  Claisen  condensation 
product,  while  only  a  trace  of  this  product  was  found  when  methyl 
formate  was  used  as  the  electrophile.  In  addition  to  the  greater 
electrophilicity  of  the  carbonyl  carbon,  o-hydrogens  are  more  acidic 
in  thiol  esters  due  to  more  effective  delocalization  of  the  negative 
charge  of  the  conjugate  base.  Hence,  thiol  esters  also  exhibit  more 
nucleophilic  activity  than  alcohol  esters.4,6 

Given,  then,  the  thiol  ester  form  of  the  substrates  and 
their  associated  nucleophilic  and  electrophilic  activity,  the 
simplest  model  for  the  condensation  in  polyketide  biosynthesis 
involves  formation  of  an  anion  and  subsequent  Claisen  type 
condensation. 8  This  proposal  was  rejected  by  Arnstadt  et  ai, 
however,  on  the  basis  of  isotopic  studies.®  In  their  studies, 
dideutero-malonyl-CoA  was  used  as  the  substrate  for  yeast  fatty  acid 
synthetase  (FAS) ,  and  no  primary  isotope  effect  on  the  rate  of  fatty 
acid  synthesis  was  observed.  Additionally,  when  they  investigated 
the  condensation  reaction  separately  from  other  FAS  reactions  by 


using  the  /9-ketoacyl -acyl -carrier -protein  synthetase  of  Esterichia 
colL  (see  below)  in  the  presence  of  tritiated  water,  they  did  not 
observe  any  incorporation  of  tritium  in  the  acetoacetyl- thiol  ester 
product.  From  these  results,  Arnstadt  concluded  free  carbanions  are 
not  involved  in  the  condensation,  and  the  formation  of  the  new 
carbon  bond  is  coupled  with  the  cleavage  of  the  carbonyl  bond  of  the 
malonyl  group. 


These  conclusions  are  based  on  a  model  in  which  enzyme 
catalyzed  reactions  take  place  in  an  analogous  manner  to  reactions 
in  solution.  Dewar  and  Storch  recently  introduced  an  alternative 
model  in  which  solvent  molecules  are  excluded  from  the  active  sites 
of  enzymes.  As  a  result,  the  enzyme  catalyzed  reactions  occur  in  a 
solvent  free  environment,  similar  to  gas  phase  reactions.10  If  this 
is  the  situation  for  the  condensation  reaction  of  polyketide 
biosynthesis,  the  participation  of  a  free  carbanion  cannot  be 
eliminated  by  the  results  of  Arnstadt  et  al,9  and  the 
decarboxylation  of  the  malonate  thiol  ester  and  condensation  with 
the  acyl  thiol  ester  may  involve  carbanions  and  may  not  be 
concerted.  Consequently,  it  was  of  interest  to  study  the 
decarboxylation  of  malonate  thiol  esters  and  malonate  alcohol 
esters,  and  their  condensation  with  acetyl  thiol  and  alcohol  esters, 
using  the  semi -empirical  quantum  mechanical  models  MNDO11  and  AMI12. 
Since  these  models  are  parameterized  using  gas  phase  experimental 
data,  results  should  reflect  what  occurs  if  the  solvent  molecules 
are  excluded  from  the  active  site  of  the  enzyme. 
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Procedure 

The  standard  closed  shell  MNDO  '  and  AMI12  procedures  as 
implemented  in  the  AMPAC14  package  of  computer  programs  were  used 
for  all  calculations.  With  the  exception  of  complexes  in  which 
separation  of  charged  groups  was  required  (see  below) ,  all 
geometries  were  optimized  by  minimizing  the  energy  with  respect  to 
all  geometrical  variables  using  the  DFP  method.15  Minima  and 
transition  states  were  characterized  whenever  possible  by  insuring 
the  applicable  Hessian  (second  derivative)  matrices  had  exactly  zero 
or  one  negative  eigenvalue,  respectively.  '  Entropies  were 
calculated1*  when  it  was  anticipated  the  entropy  of  activation  and 
the  entropy  of  reaction  would  be  significant.  The  entropy 

contribution  for  an  internal  rotation  lost  during  the  course  of  the 
reaction  was  calculated  using  the  tables  of  Pitzer  and  Gwinn.20 

Model 

The  mechanism  for  fatty  acid  synthesis  is  the  most 
thoroughly  studied  of  the  polyketide  synthesis  mechanisms.1  Since  it 
is  anticipated  the  syntheses  of  other  polyketides,  and  specifically 
the  condensation  portion  of  the  syntheses,  follow  a  pathway 
essentially  similar  to  that  of  fatty  acid  synthesis,4’6  discussion 
is  limited  to  fatty  acid  synthesis. 
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Fatty  acid  synthetases  are  found  in  two  general  forms .  In 
EsCerichia  coli  and  several  plants,  fatty  acid  synthesis  is 

catalyzed  by  at  least  six  enzymes  and  acyl  carrier  protein  (ACP) , 
nonassociated  or  loosely  associated  in  a  complex  and  separable  from 
each  other  using  conventional  techniques  for  enzyme  fractionation 
and  purification.  '  *  FAS  from  yeast  and  animals,  on  the  other 
hand,  are  multifunctional  molecular  enzyme  complexes.  '  ’  Most 
work  on  elucidating  structural  features  of  the  active  site  of  the 
condensing  activity  has  been  done  on  the  associated  type  of 
synthetase,  and  discussion  is  further  restricted  to  this  type. 


The  condensation  reaction  involves  two  active  thiol 

groups.24  One  thiol  group  belongs  to  a  cysteine  residue  localized  in 

the  condensing  activity  of  the  FAS  (hereafter  referred  to  as 

cys-SH).  The  second  thiol  group  belongs  to  the 

4' -phosphopantetheine  group  of  the  ACP  of  the  FAS  (hereafter 

referred  to  as  pant-SH) .  These  active  thiol  groups  can  approach 

within  approximately  2  A  of  each  other.2*’25  Prior  to  condensation, 

the  acetyl  transacylase  activity  of  the  FAS  transfers  an  acetyl 

"primer"  group  from  acetyl  coenzyme  A  (acetyl-CoA,  CH^OS-CoA)  to 

the  pant-SH,  from  which  it  is  transferred  again,  this  time  to  the 
-..2A.26 

cys-SH  : 
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CH  CS-CoA  +  pant-SH  -»  CH  CS-pant  +  CoA-SH 
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0  0 

*  II 

CH3CS-pant  +  cys-SH  ■+  CH3CS-cys  +  pant-SH 


The  malonyl  transacylase  activity  then  transfers  a  malonate  group 
from  malonyl  *  coenzyme  A  (malonyl -CoA,  HO^CCH^COS  -  CoA)  to  the 
pant-SH24: 


HO  CCH  CS-CoA  +  pant-SH  -*•  HO  CCH  CS-pant  +  CoA-SH 
2  2  2  2 


The  mechanisms  of  these  transacylases  are  not  known.5 
Decarboxylation  of  the  malonate  group  and  condensation  with  the 
acetyl  group,  with  inversion  of  configuration, 27  follows: 


CH3CS-cys  +  H02CCHzCS-pant  - 


0  0 
II  II 

CH  CCH  CS-pant  +  cys-SH  +  CO 


with  decarboxylation  providing  the  ultimate  driving  force  for  the 
reaction.1  The  0-ketoacyl  thiol  ester,  still  bonded  to  ACP,  then 
undergoes  reduction,  dehydration  and  reduction,  again,24  to  the 
saturated  acyl  thiol  ester.  After  transfer  of  the  acyl  thiol  ester 
to  the  cys-SH,  another  malonyl  group  is  transferred  from  malonyl -CoA 
to  the  pant-SH  and  the  condensation-reduction-dehydration-reduction 
cycle  begins  again. 
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The  « -amino  group  of  a  lysine  residue,  as  well  as  a  second 
cysteine  thiol  group,  are  also  present  in  the  active  site  of  the 
condensing  activity.28  The  e- amino  group  is  thought  to  perturb  the 
pK  of  the  active  cys-SH,  with  the  two  groups  existing  as  either  a 
hydrogen-bonded  complex  between  the  cys-SH  and  the  lys-NHz  groups,29 
an  ionic  complex  between  cys-S  and  lys-NH3  ,  ’  or  a  hybrid  of 

the  two,  in  the  absence  of  an  acyl  group  on  the  active  cysteine.  A 
mechanistic  role  for  the  lys-NH3+  group  in  the  decarboxylation  of 
HO^CCH^OS-pant  as  a  proton  donor/electron  sink  to  stabilize  the 
incipient  enolate  anion  is  also  possible.24,28  The  catalytic  role, 
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if  any,  of  the  second  cysteine  thiol  group  and  the  significance  of 
its  close  proximity  to  the  active  thiol  groups  are  unknown. 
However,  since  it  is  present  in  yeast  FAS  as  well  as  chicken  liver 
FAS,  it  also  may  play  a  role  in  the  condensation. 28 
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In  the  calculations,  a  CH3S'  group  was  used  to  model  the 
pantetheine  and  cysteine  functions  of  the  acetyl  and  malonyl  thiol 
esters.  Since  AMI  is  able  to  reproduce  hydrogen  bonds  whereas  MNDO 


is  not, 


11,12,30 


and  sulfur  has  not,  as  of  yet,  been  parameterized 


for  AMI,  corresponding  alcohol  esters  were  used  in  complexes 
involving  hydrogen  bonding.  The  «-NH2/NH3+  group  of  the  lysine 
residue  was  modeled  by  methyl  amine/ammonium  ion  rather  than 
ammonia/ammonium  ion.  This  allowed  internal  coordinates  of  the 
methyl  carbon  to  be  fixed  while  still  allowing  greater  freedom  in 
the  position  of  the  -NH3+  group  with  respect  to  the  malonyl  alcohol 
ester.  The  fixing  of  internal  coordinates  of  the  methyl  carbon  of 
methyl  ammonium  ion,  and  of  the  carboxylate  carbon  of  the  malonyl 
alcohol  ester,  was  required  in  some  calculations  to  prevent  the 
negatively-charged  carboxylate  group  and  the  positively-charged 
ammonium  group  from  collapsing.31  In  order  to  calculate  entropic 
contributions  to  the  decarboxylation,  the  required  Hessian  matrix 
had  to  be  calculated  without  constraining  any  internal  coordinates. 
For  this  purpose  formic  acid  was  used  as  the  proton  donor  for  the 
hydrogen  bond  to  the  carbonyl  group  of  the  malonyl  alcohol  ester.  A 
formic  acid  molecule  was  also  used  to  model  the  second  cys-SH.  In 
this  situation  formic  acid  was  chosen  because  a  possible  role  of  the 
second  cys-SH  may  be  in  interacting  with  the  carbonyl  oxygen  of  the 
acyl  group  bonded  to  the  active  cys-SH,  and  the  AMI  calculated 
deprotonation  enthalpy  of  formic  acid,  354  kcal/mol,30  is  close  to 
the  experimental  value  for  methane  thiol,  359  kcal/mol.32 


Table  1.  Comparison  with  experiment  of  MNDO  and  AMI  calculated  heats 
of  formation  (kcal/mol) . 
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Results 

A.  Calculation  of  Model  Compounds.  Recent  studies  have 
indicated  that,  with  few  exceptions,  AMI  and  MNDO  give  good 
estimates  of  proton  affinities  (PA)  and  deprotonation  enthalpies 
(DPE)  of  neutral  molecules.30,33  The  errors  in  PAs  and  DPEs,  as  well 
as  in  the  heats  of  formation  of  the  corresponding  cationic  acids  and 
anionic  bases,  are  similar  in  magnitude  to  those  in  the  heats  of 
formation  of  the  neutral  molecules.  Results  are  presented  in  Table 
1  for  the  calculated  heats  of  formation  of  the  model  compounds  in 
this  study.  Errors  are  generally  within  the  average  errors  for 
MNDO11  and  AMI12,  with  the  performance  of  AMI  generally  better  than 

11  33 

MNDO.  Three  previously  indicated  '  problem  areas  are  evident  for 
MNDO:  the  conjugate  acid  of  methyl  amine,  the  conjugate  base  of 

formic  acid  and  carbon  dioxide.  The  lone  significant  discrepancy 
for  AMI  is  carbon  dioxide,  which  was  also  previously  indicated.12 
Since  these  errors  are  known,  they  can  be  allowed  for  when  analyzing 
the  results  of  reaction  path  calculations.  Most  significant  in  this 
respect  are  the  errors  in  the  conjugate  bases  of  carboxylic  acids. 
MNDO  and  AMI  are  known  to  overestimate  the  DPEs  of  aliphatic 
carboxylic  acids  by  approximately  10  kcal/mol  and  6  kcal/mol, 
respectively.3033  These  values,  then,  are  used  to  estimate  errors 


for  the  malonate  derivatives  in  Table  1. 


Table  2.  Calculated  atomic  charges  for  methyl  acetate, ^methyl  thiolacetate 
and  their  conjugate  bases  resulting  from  loss  of  an  a-H  . 


H 

h^2 

/ 

H 

/3 
-  c* 

X4  -  <CH  )5 

AMI 

MNDO 

Neutral 

Anion 

Change 

Neutral  Anion 

Change 

X-0 

c1 

0.30 

0.34 

0.04 

0.35 

0.38 

0.03 

c2 

-0.22 

-0.64 

-0.42 

0.05 

-0.50 

-0.55 

o3 

-0.35 

-0.58 

-0.23 

-0.36 

-0.61 

-0.25 

o4 

-0.28 

-0.35 

-0.07 

-0.35 

-0.39 

-0.04 

cciy5 

0.21 

0.07 

-0.14 

0.23 

0.11 

-0.12 

x-s3 

c1 

0.14 

0.22 

0.08 

c2 

0.01 

-0.46 

-0.47 

o3 

-0.29 

-0.56 

-0.27 

s4 

0.03 

-0.16 

-0.19 

<CH3)S 

0.03 

-0.07 

-0.10 

The  MNDO  DPE40  for  methyl  thiolacetate  (363.6  kcal/mol)  is 
significantly  lower  than  the  MNDO,  AMI,  and  experimental  DPEs  for 
methyl  acetate  (370. 4, 40  3  6  9. 430’40  and  371. 032  kcal/mol, 

respectively).  This  is  in  agreement  with  the  trend  expected, 
considering  the  argument  concerning  more  effective  charge 
delocalization  in  the  enthiolate  anion  (see  above).  However, 
calculated  atomic  charges  (Table  2)  do  not  indicate  a  significant 
difference  in  the  amount  of  the  anionic  charge  absorbed  by  the 
carbonyl  oxygen  in  methyl  acetenthiolate  anion  as  compared  to  methyl 
acetenolate  anion.  A  larger  change  is  seen  in  the  amount  of  charge 
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absorbed  by  the  sulfur  atom  as  opposed  to  the  alcohol  oxygen  atom. 
The  calculations,  then,  model  the  difference  in  electronic  behavior 
between  the  thiol  ester  and  the  alcohol  ester  more  as  a  polarization 
effect  than  as  a  resonance  effect.  Comparison  of  the  MNDO 

calculated  heat  of  formation  of  CH^S  (-13.8  kcal/mol)  with  the 

experimental  value  (-12.2  kcal/mol32)  does  not  indicate  a 
significant  overstabilization  of  negative  charge  by  sulfur  in  MNDO 
calculations.  However,  contrary  to  conclusions  based  on  C  NMR 
results,7  the  carbonyl  carbon  is  predicted  to  carry  a  lower  positive 
charge  in  the  thiol  ester  than  in  the  alcohol  ester.  Hence, 
polarization  effects  of  the  sulfur  atom  are  probably  overestimated 
in  the  calculations,  while  resonance  effects  are  underestimated. 

B.  Decarboxylation.  Thermodynamic  results  for  various 

decarboxylation  schemes  are  given  in  Table  3 .  The  decarboxylation 

of  O^CCH^COSCH^  was  simulated  by  treating  the  C^C  •  ■  •  C^R  bon 

the  reaction  coordinate,  and  optimizing  all  other  geometric 

variables  while  incrementally  increasing  the  value  of  the  reaction 

coordinate.  The  reaction  was  predicted  to  be  endothermic  by  24.5 

kcal/mol,  with  a  forward  activation  barrier  of  29.0  kcal/mol.  The 

vibration  corresponding  to  the  negative  eigenvalue  of  the  Hessian 

matrix  of  the  transition  state  was  the  0C---CHR  stretch.  The 

2  2 

value  of  the  reaction  coordinate  was  2.41  A  in  the  transition  state. 


i 

i 
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Tsbls  3.  Thsrmodynamie  valuta  calculated  £or  decarboxylation  schemes  treating  the  0  C--CH  R 

a,b  2  2 

bond  length  as  ths  reaction  coordinate. 


* 

Reaction 

Method 

AH 

AS 

AG  AH 

AS 

AG 

rxn 

rxn 

rxn 

0 

0 

u 

■  -'i 

c 

0  CCH  CSCH  -* 

H  CCSCH  +  CO 

MNDO 

29 . 0 

3.04 

28.1  24. 5 

33.3 

14.5 

2  2  3 

2  3  2 

d 

d 

16.5 

6.3 

+ 

OH 

OH 

i 

i 

O  CCH  CSCH  - 

H  CCSCH  +  CO 

MNDO 

- 

- 

-45.8 

-55.8 

2  2  3 

2  3  2 

+ 

HNH  CH  NH  CH 

0  2 

1 

3  OH  2  3 

I 

0  CCH  COCH  - 
2  2  3 

H  CCOCH  +  CO 

2  3  2 

AMI 

10.7 

9.6  -20.9 

-30.9 

AMI 

10.  »f 

MNDO 

16. 9f 

♦ 

HNH  CH  NH  CH 

0  2 

3  OH  2  3 

i 

i 

O  CCH  CSCH  - 

H  CCSCH  +  CO 

MNDO 

13.3* 

2  2  3 

2  3  2 

HOCHO 

HOCHO 

0 

0 

M 

--'I 

O  CCH  COCH  - 

H  CCOCH  +  CO 

AMI 

- 

- 

25.4 

30.9 

16.0 

2  2  3 

2  3  2 

d 

d 

20.9 

11.6 

a 

b 

Enthalpies  and  free  energies 

are  in  kcal/mol; 

entropies 

sre  in  e.u.  T« 

-300K.  The 

barrier 

for  the  hindered  rotation  of  the  carboxylate  group  in  the  reactant  was  estimated  as  0.5 

kcal/mol,  the  barrier  given  by  Benson  for  the  CH  -CO  H  hindered  rotation.34  These  figures 

3  2 

were  calculated  using  the  experimental  AH  (CO  )  and  corrected  heats  of  formations  for 

f  2  e 

carboxylate  ions  based  on  the  errors  given  in  Table  1.  The  change  in  entropy  was  assumed  to 
be  the  same  as  that  calculated  for  decarboxylation  of  the  malonate  thiol  ester.  The  O-H 


hydrogen  bond  length  was  restricted  to  that  optimized  by  AMI  in  the  reactant  complex. 


The  large  MNDO  error  In  the  heat  of  formation  of  CO^  could 
account  for  most  of  the  heat  of  reaction  and  for  part  of  the 
activation  barrier.  However,  since  MNDO  underestimates  the 
stability  of  the  carboxylate  anion  relative  to  the  carboxylic  acid 
by  approximately  10  kcal/mol,  the  two  errors  largely  offset  each 
other.  The  result  is  that  the  reaction  is  still  predicted  to  be 
endothermic  by  a  significant  amount.  Since  these  corrections  cannot 
be  assumed  valid  for  the  heat  of  formation  of  the  transition  state, 
no  corrected  value  can  be  estimated  for  the  activation  barrier.  In 
any  case,  it  is  clear  the  activation  barrier  is  formidable. 

With  the  gain  of  three  translational  and  three  rotational 
degrees  of  freedom  at  the  expense  of  five  vibrations  and  a  hindered 
internal  rotation,  the  decarboxylation  will  be  favored  entropically . 
The  increase  in  entropy  is  not  large  enough,  however,  to  offset  the 
unfavorable  enthalpy  changes.  It  is  reasonable  to  conclude,  then, 
that  there  must  be  other  interactions  with  the  malonate  thiol  ester 
to  prompt  decarboxylation. 

As  discussed  above,  the  carbonyl  group  beta  to  the 
carboxylate  group  is  expected  to  activate  the  malonate  thiol  ester 
toward  decarboxylation  by  allowing  the  delocalization  of  the 
negative  charge  of  the  product  enthiolate  anion.  Jencks  has  pointed 
out  that  the  carbonyl  group  is  an  even  more  active  electron  sink 
when  it  is  protonated. 4 1  As  expected  on  this  basis,  when  a  proton 
was  added  to  the  carbonyl  group  and  the  geometry  optimized,  CO ■  was 
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lost  without  activation.  The  decarboxylation  was  exothermic  by  45.8 
kcal/mol  and  exergonic  by  55.8  kcal/mol. 


A  logical  source  of  this  proton  would  be  an  ammonium  ion 

derived  from  the  e-NH2  group  of  the  lysine  residue  in  the  active 

site.  To  simulate  it,  methyl  ammonium  ion  was  positioned  in  the 

plane  of  the  ester  linkage  of  the  malonate  alcohol  ester,  allowing 

hydrogen  bonding  between  the  ammonium  ion  and  the  ester  function. 

Bond  angles  and  dihedral  angles  were  frozen  to  prevent  the  ammonium 

ion  and  the  carboxylate  group  from  collapsing.  The  decarboxylation 

was  now  found  to  be  exothermic  and  exergonic  (-20.9  and  -30.9 

kcal/mol,  respectively).  An  apparent  transition  state17  at  an 

0  C---CHR  distance  of  2.06  A  defined  an  activation  barrier  of  10.7 
2  2 

kcal/mol.  A  proton  was  transferred  from  the  methyl  ammonium  ion  to 
the  carbonyl  group  of  the  ester  at  a  C^C-.-CH^R  distance  of  3.9  A. 
Interestingly,  however,  when  the  errors  for  carboxylate  anions  and 
CO  were  taken  into  account,  the  enthalpy  of  reaction  for  the 
decarboxylation  was  approximately  zero  prior  to  proton  transfer  (not 
shown  in  Table  3).  Furthermore,  when  the  bond  between  the  proton  on 
the  methyl  ammonium  ion  and  the  carbonyl  oxygen  was  treated  as  the 
reaction  coordinate,  the  proton  was  transferred  to  the  carbonyl 
(AH**-  12.2  kcal/mol;  AH  -  8.2  kcal/mol)  without  resulting  in 

rxn 

decarboxylation  (see  below) . 


.1 AJ  j.  I  «  W5W«  *. 


TO  W  R IJ  VJ  VC  W5IW7!  v.  r.  \pjvt*.wjv '.  v:  r.  r; 


Page  104 

This  model  was  oversimplified  in  three  respects.  First,  it 
ignored  possible  electrostatic  interactions  with  the  electrophile 
(thiol  ester)  with  which  the  nascent  carbanion  reacts  in  the  next 
step.  This  will  be  addressed  later  when  discussing  the  possibility 
of  a  concerted  decarboxylation- condensation  process.  Secondly.it 
replaced  oxygen  for  sulfur  in  the  ester  linkage.  Thirdly,  it 
restricted  some  of  the  internal  coordinates  of  the  malonate  ester 
ammonium  ion  complex. 

To  evaluate  the  effect  of  substituting  oxygen  for  sulfur  in 
the  ester  linkage,  the  hydrogen  bond  length  was  fixed  at  the  AMI 
optimized  value  during  decarboxylation.  Decarboxylation  was  then 
simulated  using  AMI  and  MNDO  with  the  alcohol  ester  complex,  and 
using  MNDO  with  the  thiol  ester  complex.  Three  contributions  to 
changes  in  the  reaction  profile  were  taken  into  account  when 
evaluating  the  results:  (1)  the  effect  of  restricting  the  hydrogen 
bond  length,  (2)  the  effect  of  using  MNDO  instead  of  AMI,  and  (3) 
the  effect  of  having  oxygen  rather  than  sulfur  in  the  ester  linkage. 
Restricting  the  hydrogen  bond  length  increased  the  barrier  only  from 
10.7  to  10.9  kcal/mol.  Using  MNDO  rather  than  AMI  increased  the 
barrier  another  6.0  kcal/mol.  This  is  consistent  with  previous 

comparisons  of  AMI  and  MNDO  calculated  activation  barriers . 12,42 
Additionally,  the  MNDO  barrier  for  the  alcohol  ester  was  3.6 
kcal/mol  larger  than  that  for  the  thiol  ester.  Using  these 
corrections,  the  barrier  to  decarboxylation  of  the  malonate  thiol 
ester  with  an  optimized  hydrogen  bond  length  can  be  estimated  as  7.1 


7.T.7: 


•>.v; 
•  y.v. 

§  v-'-. 


W.vv 


•>.y 


etZxdbtm.*?, A  «* -  t-  a'a*^.V.V.  •  / 


kcal/mol. 


Next  considered  was  the  effect  of  the  restrictions  on  bond 
angles  and  dihedral  angles  of  the  methyl  ammonium  ion  and  the 
carboxylate  group,  which  were  needed  to  prevent  the  oppositely 
charged  groups  from  reacting  with  one  another.  To  evaluate  this, 
formic  acid  was  used  to  form  a  hydrogen  bond  in  place  of  methyl 
ammonium  ion,  and  the  geometry  was  fully  optimized.  Reactants  and 
products  of  the  decarboxylation  were  characterized  as  true  minima. 
While  the  decarboxylation  was  endothermic,  the  enthalpy  increased 
smoothly  throughout  the  reaction,  with  no  transition  state  between 
reactant  and  products.  It  is  therefore  possible  that  at  least  part 
of  the  apparent  activation  barrier  in  the  decarboxylation  of  the 
malonate  ester  -  ammonium  ion  complex,  and  certainly  the 
corresponding  thiol  ester  complex,  is  an  artifact  due  to 
restrictions  placed  on  the  geometry  in  the  calculations. 

The  reaction  path  calculated  for  the  decarboxylation  of  the 
malonate  ester  -  ammonium  ion  complex  ultimately  involved  proton 
transfer.  When  the  proton  transfer  from  the  methyl  ammonium  ion  to 
the  malonate  ester  was  treated  as  the  reaction  coordinate,  there  was 
an  activation  barrier  of  12.2  kcal/mol.  The  question  remains,  then, 
as  to  whether  the  apparent  barrier  to  decarboxylation  might  be 
partly  due  to  the  barrier  for  proton  transfer.  There  is  evidence, 
however,  that  this  is  not  the  case.  First,  the  barrier  to 
decarboxylation  occurs  much  earlier  along  the  reaction  profile  than 


the  proton  transfer.  In  the  apparent  transition  state  for  the 
decarboxylation,  the  H-N  bond  for  the  proton  ultimately  transferred 
to  the  carbonyl  oxygen  of  the  ester  has  stretched  only  0.01  A  from 
what  it  was  in  the  reactant  complex.  Additionally,  restricting  the 
length  of  the  0---H  hydrogen  bond  increased  the  barrier  to 
decarboxylation  by  only  0.2  kcal/mol  (see  above).  Finally,  proton 
transfer  alone  did  not  cause  decarboxylation.  These  results 
indicate  that  the  barrier  to  proton  transfer  plays  little  or  no  role 
in  the  decarboxylation. 43 

This  raises  yet  another  question:  Is  proton  transfer 
required  at  all,  or  is  the  electrostatic  interaction  between  the 
malonate  thiol  ester  and  the  ammonium  ion  alone  sufficient  to 
initiate  decarboxylation?  When  the  ammonium  ion  of  the  malonate 
ester  -  methyl  ammonium  ion  -  methyl  acetate  complex  was  replaced  by 
a  pure  ionic  charge  1.85  A  from  the  carbonyl  oxygen  of  the  malonate 
(the  same  distance  as  NH---0-C  just  prior  to  proton  transfer 

malonate 

during  decarboxylation) ,  the  enthalpy  of  activation  for 
decarboxylation  was  7.2  kcal/mol.  When  the  separation  was  increased 
to  2.25  A,  a  distance  large  enough  to  inhibit  proton  transfer,  the 
barrier  to  decarboxylation  increased  only  to  12.0  kcal/mol.  The 
correction  for  using  the  alcohol  ester  rather  than  the  thiol  ester 
would  decrease  this  barrier  to  less  than  10  kcal/mol.  It  is 
feasible,  therefore,  the  charge  of  the  ammonium  ion  alone  would 
prompt  decarboxylation  in  the  enzyme. 
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As  noted  above,  it  has  been  claimed  that  free  carbanions 

cannot  be  involved  in  the  enzyme  reaction  because  of  the  lack  of 

tritium  incorporation  when  it  is  carried  out  in  tritiated  water.9 

Such  exchange  will,  however,  be  impossible  if  water  is  excluded  when 

the  malonate  is  absorbed  into  the  active  site.  To  check  the  effect 

of  water,  five  water  molecules  were  added  to  the  model  system. 

Whereas  the  protonated  alcohol  ester  and  the  protonated  thiol  ester 

both  lost  CO  without  activation  in  the  absence  of  water,  the 
2 

"hydrated"  protonated  malonate  ester  did  not  dissociate.  When  the 
O^C-'-CH^R  bond  length  was  treated  as  the  reaction  coordinate,  a  5.7 
kcal/mol  barrier  to  decarboxylation  was  predicted  by  AMI.  Given 
that  the  barrier  would  be  much  larger  for  a  mechanism  involving 
interaction  of  the  0-carbonyl  group  with  an  ammonium  ion,  and  that 
water  molecules  could  also  weaken  this  interaction,  these  results 
support  the  suggestion  that  a  crucial  factor  in  enzyme  reactions  is 
the  exclusion  of  water  from  the  reacting  system. 


The  decarboxylation  of  the  malonate  thiol  ester  in  a  solvent 
free  environment  thus  appears  to  be  strongly  exothermic  and 
exergonic  when  an  ammonium  ion  is  near  enough  for  proton  transfer  to 
occur  during  the  reaction.  However,  even  without  proton  transfer, 
the  enthalpy  of  reaction  is  approximately  zero.  Since  the  entropy 
of  decarboxylation  is  negative,  the  free  energy  change  for  the 
reaction  without  proton  transfer  will  be  negative.  Furthermore, 
there  is  a  low  enthalpic  barrier  to  decarboxylation  under  these 
conditions  (<10  kcal/mol) . 


»S 
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C.  Condensation.  The  second  step,  involving  Claisen 
condensation  of  the  anion  and  ester,  was  studied  by  treating  the 
distance  between  the  methene  carbon  of  the  enolate  or  enthiolate 
anion  and  the  carbonyl  carbon  of  the  corresponding  alcohol  or  thiol 
ester  as  the  reaction  coordinate .  As  the  value  of  the  reaction 
coordinate  was  incrementally  decreased,  all  other  internal 
coordinates  were  fully  optimized,  except  where  a  few  internal 
coordinates  were  frozen  to  prevent  the  collapsing  of  oppositely 
charged  groups.31  Table  4  lists  the  enthalpy  of  activation  and 
enthalpy  of  reaction  results .  Entropy  changes  will  not  play  a 
significant  role  in  the  corresponding  enzymatic  reaction  since  all 
intermediates  remain  bound  to  the  enzyme.  Consequently,  entropies 
were  not  calculated. 

Condensation  of  the  enthiolate  anion  with  methyl 
thiolacetate  was  predicted  by  MNDO  to  be  2.6  kcal/mol  exothermic, 
with  an  activation  barrier  of  15.3  kcal/mol.  Given  the  propensity 
of  MNDO  to  overestimate  activation  barriers,12  the  condensation  of 
the  oxygen  analogs  was  simulated  using  MNDO  and  AMI.  MNDO  predicted 
essentially  the  same  barrier,  14.5  kcal/mol.  AMI,  on  the  otherhand, 
predicted  the  exothermic  formation  of  a  charge -dipole  complex  at  4.5 
A  (AH  -  -7.4  kcal/mol),  which  was  characterized  as  a 

compilation 

true  minimum.  From  the  charge-dipole  complex  there  was  a  5.8 
kcal/mol  barrier  to  formation  of  the  tetrahedral  adduct.  The  lower 
barrier  calculated  by  AMI  as  compared  to  MNDO  is  consistent  with 
previous  results  for  the  formation  of  tetrahedral  complexes.42  Since 
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the  positive  charge  on  the  carbonyl  carbon  of  the  thiol  ester  is 
underestimated  by  MNDO,  the  barrier  to  condensation  for  the  sulfur 
compounds  should  be  at  least  as  low.  Additionally,  with  the 
enthiolate  anion  formed  while  bonded  in  the  active  site  of  the 
enzyme  in  a  complex  tight  enough  to  exclude  water  molecules,  it  may 
be  formed  closer  to  the  methyl  thiolester  than  the  4.5  A  of  the 
methyl  acetate  -  methyl  acetenolate  complex.  In  this  case,  the 
barrier  to  condensation  of  the  enthiolate  anion  and  thiol  ester  in 
the  enzyme  would  be  less  than  5 . 8  kcal/mol  from  the  point  on  the 
reaction  coordinate  at  which  the  enthiolate  anion  is  formed. 

As  modelled  by  AMI,  however,  the  decarboxylation  of  the 
malonate  ester  ultimately  involves  proton  transfer  from  the  methyl 
ammonium  ion  to  the  carbonyl  oxygen  of  the  enolate  anion.  The 
condensation,  in  this  case,  would  involve  an  enol  rather  than  an 
enolate  anion.  Without  proton  transfer  simultaneous  with  the 
condensation,  a  zwitterionic  product,  1,  would  be  formed: 
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Tablt  A.  Thermodynamic  values  calculated  for  condensation  schemes 

treating  the  RCH  •••  OO  distance  as  the  reaction  coordinate. 
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kcal/mol.  Values  based  on  the  heats  of  formation  of  the  isolated 
c 

reactants.  Values  based  on  the  heat  of  formation  of  the  charge-dipole 

d 

complex  between  the  reactants  (see  text).  A  transition  state  was  not 

e 

found  (see  text).  This  value  differs  by  a  few  kcal/mol  depending  on  the 
finel  distance  of  the  amine  from  the  hemiacetal,  in  this  case  2.8  A  from 
the  hydroxyl  hydrogen. 
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In  Che  enzyme,  the  proton  transfer  could  be  envisaged  in 
either  of  two  ways:  (1)  transfer  directly  from  the  carbonyl  oxygen 
of  the  enol  moiety  to  the  carbonyl  oxygen  of  the  acetate  moiety  via 
a  six  membered  cyclic  transition  state,  or  (2)  proton  transfer  from 
the  enol  moiety  to  the  amine  and  proton  transfer  from  the  amine  to 
the  methyl  acetate  moiety.  The  latter  picture  is  essentially  a 
proton  relay  model,  and  is  analogous  to  the  "charge  transfer 
mechanism”  of  chymo trypsin. 46  Both  of  these  processes  would  give  the 
same  net  result,  a  hemiacetal  and  methyl  amine,  and  are  exothermic. 
A  grid  search  for  the  direct  proton  transfer  during  condensation 
indicated  the  enthalpy  of  activation  would  be  of  the  order  of  30 
kcal/mol.  Given  the  size  of  this  barrier  and  the  proximity  of  the 
amine  group  to  the  condensing  moieties  in  the  enzyme,  as  well  as  its 
involvement  to  this  point  in  the  reaction,  direct  proton  transfer  is 
not  likely. 

To  obtain  a  geometry  for  the  complex  involved  in  the 
proton- relay  model,  the  amine  group  was  positioned  to  interact  with 
the  carbonyl  groups  of  the  raalonate  ester  and  methyl  acetate. 
Constraining  only  the  bond  angle  and  dihedral  angle  of  the 
carboxylate  carbon  to  keep  the  carboxylate  group  approximately 
perpendicular  to  the  plane  of  the  ester  linkage,31  the  complex 
optimized  at  a  separation  of  3.68  A  between  the  two  carbons  which 
would  form  the  new  C-C  bond  upon  condensation  (N-H---0-C 

malonate 

2.12  A).  When  the  carboxylate  group  was  removed  and  the  geometry 
was  optimized  with  the  carbon-carbon  distance  fixed  at  3.68  A,  a 
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proton  transferred  from  the  amine  to  the  carbonyl  oxygen  of  the 
enolate  to  form  the  enol.  The  carbon- carbon  distance  was  then 
treated  as  the  reaction  coordinate  and  incrementally  decreased  while 
optimizing  all  other  internal  coordinates.  The  enthalpy  increased 
steadily  until  the  apparent  intermediate  3  was  obtained  30.1 
kcal/mol  above  the  initial  complex.  A  force  constant  analysis 
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indicated  3  was  not  a  true  minimum.  Instead,  there  was  one 
negative  eigenvalue  corresponding  to  movement  of  the  proton  from  the 
enol  moiety  to  the  amine.  The  resulting  complex,  4,  was  also  not  a 
minimum,  optimizing  instead  to  5.  The  hemiacetal,  5,  then,  was 
characterized  as  the  only  true  minimum  resulting  from  formation  of 
the  new  carbon-carbon  bond  between  the  enol  and  methyl  acetate. 
Hence,  the  "true"  reaction  coordinate  involved  not  only  the 
formation  of  the  carbon-carbon  bond,  but  also  a  double  proton 
transfer.  However,  all  efforts  to  find  the  transition  state  along 
this  reaction  coordinate  failed  (see  below) . 
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D.  Concerted  Decarboxylation-Condensation.  To  detect  a 
coupling  between  these  reactions,  grid  searches  were  calculated 
using  various  geometries  of  the  malonate  thiol  ester  -  methyl  thiol 
acetate  complex  and  the  oxygen  analog,  as  well  as  the  malonate  ester 
-  methyl  ammonium  ion  -  methyl  acetate  complex.  The  O^C-'-C^R  bond 
length  and  the  separation  between  the  carbons  forming  the  new  bond 
were  treated  as  the  two  reaction  coordinates.  In  all  cases, 
decarboxylation  occurred  essentially  independently  of  condensation. 
Any  path  involving  significant  breaking  of  the  bond  to  the 
carboxylate  group  and  simultaneous  formation  of  the  new 
carbon-carbon  bond  involved  an  activation  barrier  of  at  least  50 
kcal/mol . 

It  seemed  possible  that  the  second  cysteine  in  the  active 
site  of  the  condensing  enzyme  might  participate  by  interacting  with 
the  carbonyl  group  of  the  acetate,  increasing  its  electrophilicity . 
However ,  when  the  calculation  was  repeated  with  an  added  molecule  of 
formic  acid  to  model  this  interaction,  results  were  not 
significantly  different. 

The  presence  of  methyl  acetate  did  facilitate 
decarboxylation  of  the  malonate  ester.  This  effect  was  manifested 
in  two  ways.  When  the  O^C'-'CH^R  distance  was  treated  as  the 
reaction  coordinate  in  the  malonate  ester  -  methyl  ammonium  ion  - 
methyl  acetate  complex  described  above  (C---C  distance  fixed  at  3.68 
A),  the  enthalpy  of  activation  was  only  5.2  kcal/mol,  as  opposed  to 


10.7  kcal/mol  for  the  malonate  ester  -  methyl  ammonium  ion  complex. 
Secondly,  although  there  was  again  no  proton  transfer  from  the 
ammonium  ion  to  the  oxygen  of  the  incipient  enolate  anion  until  well 
after  the  transition  state  for  decarboxylation  was  passed,  when  the 
NH  •  •  • o— C  distance  was  treated  as  the  reaction  coordinate, 

malonate 

proton  transfer  did  result  in  decarboxylation.  It  was  possible  this 
effect  resulted  from  the  different  relative  position  of  the  ammonium 
ion  in  the  malonate  ester  -  methyl  ammonium  ion  -  methyl  acetate 
complex  as  opposed  to  in  the  malonate  ester  -  methyl  ammonium  ion 
complex,  rather  than  from  the  presence  of  methyl  acetate.  To  test 
this  possibility,  the  bond  angle  of  the  nitrogen  atom  of  the 
ammonium  ion  relative  to  the  malonate  ester  was  frozen  and  the 
methyl  acetate  was  removed  from  the  calculation.  The  dihedral  angle 
of  the  nitrogen  atom  did  not  change  significantly  and  was  allowed  to 
optimize.  This  time  decarboxylation  did  not  occur  upon  proton 
transfer.  Hence,  the  facilitation  of  decarboxylation  was  the  result 
of  the  presence  of  methyl  acetate,  not  the  change  in  the  relative 
position  of  the  ammonium  ion.  Assuming  the  same  correction  as 
before  for  using  the  alcohol  ester  instead  of  the  thiol  ester,  and 
the  same  value  for  the  entropy  of  activation  for  decarboxylation, 
the  free  energy  of  activation  for  decarboxylation  can  be  estimated 
as  less  than  1  kcal/mol.  This  leaves  little  doubt  decarboxylation 
would  occur  spontaneously  in  the  malonate  thiol  ester  -  ammonium  ion 
-  thiol  acetate  complex  of  the  condensing  enzyme. 


Discussion 


The  concerted  decarboxylation- condensation  model  of  the 
condensation  reaction  in  polyketide  biosynthesis  was  proposed  by 
Amstadt  et  al  due  to  the  lack  of  a  primary  deuterium  isotope  effect 
when  dideutero-malonyl-CoA  was  used  as  the  substrate  for  fatty  acid 
synthetase,  and  the  lack  of  tritium  incorporation  when  the 
condensation  was  carried  out  in  tritiated  water.9  These  experiments 
were,  however,  of  what  might  be  termed  "one-way"  type10*  in  the 
sense  that  while  one  of  the  two  possible  outcomes  would  have  been 
significant,  the  other  is  not.  Thus  if  isotopic  exchange  had  been 
observed,  this  would  have  provided  strong  evidence  for  participation 
by  ionic  intermediates.  Failure  to  observe  it  could,  however,  be 
due  to  the  factor  noted  above,  i.e.  that  the  special  characteristic 
of  enzyme  reactions  probably  arise  from  the  fact  that  absorption  of 
the  substrate  into  an  active  site  involves  displacement  of  solvent 
(water)  from  between  them.  Here,  even  if  anions  are  involved,  they 
clearly  cannot  undergo  isotopic  exchange  with  solvent  if  there  is  no 
solvent  present. 

As  noted  above,  our  calculations  indicate  that  the 
decarboxylation  and  condensation  cannot  take  place  as  a  simple 
synchronous  process  because  the  calculated  activation  energy  was 
much  too  large.  This  result  could  have  been  anticipated  because  the 
S£2  reaction,  being  an  autoactivated  process  due  to  the  steric 
difficulties  associated  with  forming  a  pentavalent  adduct  of 
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carbon,  cannot  act  as  a  driving  force  for  formation  of  the 
reactants  in  it. 

On  the  other  hand,  the  calculations  do  imply  that  the 
condensation  and  decarboxylation  are,  at  least,  weakly  concerted. 
The  precise  geometry  of  the  malonate  thiol  ester,  lysine  residue  and 
thiol  acetate  in  the  active  site  of  the  condensing  enzyme  is  not 
known.  It  is  possible  the  methene  carbon  of  the  malonate  thiol 
ester  is  close  enough  to  the  carbonyl  carbon  of  the  thiol  acetate  at 
the  time  decarboxylation  occurs  that  there  is  even  greater  coupling 
between  the  two  reactions. 

A  model  for  the  decarboxylation- condensation  reaction  of 
polyketide  synthesis  consistent  with  experimental  evidence,  then, 
involves  the  initial  transfer  of  an  acetyl  group  from  acetyl -CoA  to 
the  active  cysteine.  The  aqueous  pKs  of  the  thiol  group  of  cysteine 
and  the  e-NH*  group  of  lysine  are  8.18  and  10.53,  respectively, 50 
while  the  gas  phase  DPEs  of  alkyl  thiols  are  on  the  order  of  350 
kcal/mol  and  the  PAs  of  alkyl  amines  are  approximately  220 
kcal/mol  With  the  close  proximity  of  the  amino  group  perturbing 

the  acidity  of  the  thiol  group,  the  lysine  and  cysteine  residues 
will  exist  as  the  ammonium  salt  in  the  presence  of  water.  The 
energy  required  to  extrude  the  water  during  the  preliminary  stages 
of  the  acylation  process,  then,  is  provided  by  the  exothermicity  of 
the  reverse  proton  transfer  in  the  now  solvent  free  environment. 
The  association  of  the  amino  and  thiol  groups  remains  close, 
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however,  and  the  proton  Is  transferred  back  to  the  amino  group 
during  the  actual  acylation  of  the  thiol  group. 

A  malonyl  group  is  then  transferred  from  malonyl-CoA  to  the 
pant-SH.  Whether  the  malonyl  group  is  deprotonated  after  binding  in 
the  malonyl  transferase  or  not  until  after  transferring  to  the 
pant*SH  is  not  known.  However,  the  requirement  for  a  solvent  free 
environment  for  the  decarboxylation  dictates  the  malonyl -CoA  be 
protonated  when  absorbed  from  solution.10*’15  In  the  presence  of  the 
acylated  cysteine  residue  and  the  protonated  lysine  residue,  the 
deprotonated  malonyl  group  decarboxylates  and  condenses  with  the 
acetyl  carbonyl  carbon.  The  resulting  hemithiolacetal ,  6, 
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then  rearranges  to  the  0-ketoacyl  thiolester  and  the  initial 
cys-SH- • • lys-NHz  hydrogen  bonded  complex. 

The  inability  to  find  a  transition  state  for  the  proton 
relay  condensation  mechanism  prevents  a  conclusive  evaluation  of 
this  model,  but  does  not  eliminate  this  model  as  a  possibility.  The 
difficulty  in  finding  the  transition  state  for  the  proton  relay 
mechanism  could  be  the  result  of  one  or  a  combination  of  several 
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factors,  both  computational  and  chemical.  Since  several  elementary 
processes  are  involved  (i.e.  formation  of  a  carbon- carbon  bond  and 
two  hydrogen  transfers) ,  an  extensive  grid  search  involving  a 
separate  coordinate  for  each  process  may  be  required  to  locate  the 
transition  state.  Such  a  three  dimensional  grid  search  involving 
the  model  system  in  this  study  is  not  feasible. 


The  problem  of  locating  a  transition  state  may  also  have  a 
chemical  origin.  As  mentioned  above,  the  actual  geometries  of  the 
substrates  in  the  enzyme  are  not  known,  and  the  possibilities  are 
virtually  endless.  Of  course,  the  location  of  a  transition  state 
will  ultimately  depend  on  selection  of  the  proper  geometries.  While 
a  large  number  of  possibilities  were  tried,  it  is  possible  the 
"correct"  one  was  not.  Especially  significant  is  the  separation  of 
the  condensing  carbons  upon  decarboxylation  of  the  malonate  thiol 
ester.  If  they  are  close  enough,  decarboxylation  and  condensation 
may  be  coupled  enough  to  prevent  any  proton  transfer.  Additionally, 
although  AMI  typically  calculates  activation  barriers  more 
accurately  than  MNDO,  it  is  known  to  seriously  overestimate  barriers 
to  some  proton  transfers.12'*9  In  this  case,  that  difficulty  would 
be  compounded  by  the  reaction  coordinate  involving  two  proton 
transfers.  Finally,  at  some  point  along  the  reaction  coordinate  a 
double-welled  potential  may  develop,  and  proton  tunnelling  may  play 
a  significant  role.  In  any  of  these  cases,  the  transition  state  as 
visualized  above  may  be  extremely  high  in  energy  on  the  AMI  PE 
surface,  and  very  difficult  to  locate. 


Based  on  the  results  above,  however,  there  is  a  simpler 
model  which  eliminates  the  above  difficulties.  As  modelltd,  proton 
transfer  does  not  occur  until  well  after  the  transition  state  for 
decarboxylation  is  passed,  and  does  not  appear  to  play  a  significant 
role  in  the  decarboxylation,  aside  from  making  the  overall  reaction 
more  exothermic.  While  it  is  known  the  lys-NH2  is  close  enough  to 
the  cys-SH  and  pant-SH  to  interact  with  them,24  it  is  not  known  how 
close  it  is  to  the  carbonyl  oxygens  of  bonded  acyl  groups.  As  such, 
the  effect  of  the  ammonium  group  in  the  enzyme  may  be  purely  an 
electrostatic  effect.  The  protonated  amine  group  may  be  far  enough 
away  from  the  carbonyl  group  of  the  malonate  that  its  positive 
charge  initiates  decarboxylation  without  actually  transferring  a 
proton  in  the  process .  Hence ,  the  decarboxylation  and  condensation 
may  take  place  without  actual  proton  transfers  except  from  the 
cys-SH  to  the  lys-NHz  and  back.  In  this  case,  condensation  should 
take  place  with  only  a  small  activation  barrier,  as  with  the  enolate 
and  acetate  alone. 

This  further  illustrates  the  difference  between  enzymatic 
reactions  and  reactions  in  aqueous  solution.  In  the  solvent  free 
environment  of  the  active  sites  of  enzymes,  any  charge  present  will 
have  a  far-reaching  effect  on  reactions.  In  solution,  this  effect 
would  be  severely  attenuated  by  the  dielectric  properties  of  water. 


Both  of  the  above  models  are  consistent  with  the  inhibition 
of  the  condensing  activity  by  iodoacetamide .  Iodoacetamide  is  known 
to  acylate  the  cys-SH:50 

lys-NH-.  lys-NH  + 

2  > 

ICH  CONH  +  cys-S  -  cys-SCH  CONH  +  I- 

2  2  J  J  2  2 

Whereas  unmodified  fatty  acid  synthetase  in  the  absence  of 
acetyl-CoA  and  NADPH  decarboxylates  malonyl  residues  very  slowly, 
the  enzyme  inhibited  by  iodoacetamide  decarboxylates  malonyl -CoA 
much  more  rapidly,  while  at  the  same  time  preventing  condensation. 50 
Kresze  et  al  attributed  this  effect  to  the  conf irmational  change  in 
the  enzyme  resulting  from  the  binding  of  the  iodoacetamide.  While 
this  may  be  part  of  the  explanation,  the  binding  of  the 

iodoacetamide  also  results  in  the  formation  of  the  free  ammonium  ion 
required  for  decarboxylation.  The  carbonyl  carbon  of  the 

iodoacetamide  is  less  electrophilic  than  that  of  an  acyl  thiol 

ester,  and  the  tetrahedral  condensation  product  probably  does  not 

form  to  a  significant  extent.  As  a  result,  the  enol  or  enolate 
eventually  picks  up  a  proton  and  is  eliminated  as  acetyl-CoA, 
opening  the  way  for  the  decarboxylation  of  another  malonate  residue. 
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ordering  of  the  calculated  heats  of  reaction  were  nearly  identical 
when  C-l-NP  and  TNT  were  the  reactants,  providing  hope  that  the 
energetic  ordering  of  the  transition  states  of  the  model  system 
would  also  parallel  those  of  TNT.1  While  experimental  evidence  for 
TNT  was  most  consistent  with  intermolecular  hydrogen  transfer,1  the 
calculated  activation  energies  for  all  C-l-NP  reactions  except  NC>2 
dissociation  were  reasonably  close  to  the  observed  activation  energy 
for  TNT  thermolysis.  The  calculated  R_N02  bond  dissociation 
enthalpy  (BDE) ,  65  kcal/mol,4  was  15-30  kcal/mol  larger  than  the 
MINDO/3  calculated  activation  enthalpies  for  the  other 
intramolecular  processes. 

Dewar  et  al  similarly  found  simple  R-N02  bond  scission  to  be 
unfavorable  in  the  thermolysis  of  nitromethane . 5  Again  using 
MINDO/3,  they  showed  that  a  nitro  to  nitrite  rearrangement,  followed 
by  dissociation  of  NO,  was  more  consistent  with  gas  phase  kinetic 

data  than  was  the  direct  dissociation  of  NO  . 

2 

These  findings  are  in  direct  contrast  to  gas  phase  pyrolysis 
results  for  substituted  nitrobenzenes . 9  Under  conditions  in  which 
surface  catalyzed  reactions  were  minimized,  Gonzalez  et  al 
determined  the  major  process  involved  scission  of  the  Ar-N02  bond, 
although  the  nitro  to  nitrite  rearrangement  and  subsequent  NO 
dissociation  also  occurred  to  a  lesser  extent.  This  was  true  for 
all  compounds,  even  those  with  a  methyl  group  ortho  to  the  departing 
nitro  group,  and  despite  the  fact  the  R-N02  BDE  derived  by  Gonzalez 
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(70  kcal/mol)  was  even  larger  than  that  calculated  by  Turner  and 
Davis.1  Indeed,  the  largest  Arrhenius  preexponential  factor  and 
smallest  Arrhenius  activation  energy  derived  from  the  nitrobenzene 
results  were  for  an  ortho  methyl  derivative,  and  this  effect  was 
attributed  to  the  steric  effects  of  the  methyl  group.  At  the  same 
time,  there  was  no  evidence  of  direct  participation  of  the  methyl 
group  in  the  decomposition. 


The  dissociation  of  NO  and  the  nitro  to  nitrite 

2 

rearrangement  with  subsequent  NO  dissociation  are  known  to  also  play 
a  role  in  the  decomposition  of  molecular  ions  of  nitro  compounds. 
Both  decomposition  channels,  for  example,  have  recently  been  shown 
to  be  active  in  the  fragmentation  of  nitromethane7  and  nitrobenzene8 
molecular  ions. 


In  the  mass  spectral  fragmentation  of  the  C-l-NP  molecular 
ion  (C-l-NP*),  these  reactions  would  compete  with  fragmentation 
involving  direct  participation  of  the  methyl  group.  Given  the  above 
theoretical  and  experimental  results  for  neutral  nitro  compounds,  it 
is  of  interest  to  determine  if  methyl  group  participation  in  the 
fragmentation  of  C-l-NP*  does  occur,  and  to  compare  the  calculated 
enthalpic  requirements  of  the  competing  processes. 


SOCK 
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Tabla  1.  Ralativa  abundances  in  the  70  eV  El  mass 

a,b 

spectra  of  cis-  and  trans- 1-nitropropene. 
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cis 

trans 
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cis 

trans 

15 

6.7 

3.6 

41 

53.9 

40.9 

26 

7.8 

3.9 

42 

43.5 

4.3 

29 

6.4 

3.9 

43 

11.0 

24.7 

30 

17.1 

12.2 

44 

13.6 

33.1 

31 

6.4 

0.2 

46 

5.5 

4.6 

37 

7.2 

5.1 

70 

21.2 

0.7 

38 

11.6 

9.4 

74 

6.4 

1.2 

39 

100.0 

89.7 

86 

11.3 

1.9 

40 

9.0 

6.2 

87 

21.4 

100.0 

a  b 

Spectra  from  ref. 9  Peaks  are  listed  only  if  the 

relative  abundance  for  at  least  one  of  the  two  isomers 
c 

exceeds  5.0.  Mass  to  charge  ratio. 

Additionally,  there  is  significant  fragmentation  of  C^O  +  N0+  (m/z 
30)  and  C3H3+  +  N02  (m/z  41)  in  both  isomers.  Loss  of  NO  and  N02+ 
(m/z  57  and  m/z  46,  respectively)  appears  to  be  relatively  minor, 
assuming  negligible  secondary  fragmentation. 


The  spectra  clearly  have  differences,  however.  This 
indicates  cis-trans  isomerization  does  not  occur  to  a  significant 
extent,  a  somewhat  surprising  result  considering  AMI  predicts  the 
barrier  to  rotation  around  the  C1-Cz  bond  to  be  only  11  kcal/mol.10 
The  spectra  also  show  more  extensive  fragmentation  of  C-l-NP*  than 
the  T-l-NP  molecular  ion  (T-l-NP*).  This  is  evident  from  the 


difference  in  the  relative  abundances  of  the  molecular  ions,  and 
from  the  larger  number  of  significant  peaks  (relative  abundance 
>5.0)  in  the  C-l-NP  spectrum.  Of  particular  note  is  the  presence  of 
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an  m/z  70  peak  in  the  C-l-NP  spectrum,  presumably  from  the  loss  of  a 
hydroxyl  radical  from  C-l-NP*,  and  its  absence  from  the  T-l-NP 
spectrum.  The  implication  of  these  results  is  that  there  is  direct 
methyl  group  participation  in  the  fragmentation  of  an  m/z  87  ion. 

As  will  be  seen,  the  m/z  86  peak  is  also  significant  in  this  regard. 

Calculation  Procedure 

The  MNDO,  MIND0/3 ,  and  AMI11  semiempirical  quantum 
mechanical  models  have  been  shown  to  give  results  comparable  to  ab 
initio  models  using  split  basis  sets  in  a  fraction  of  time  required 
for  the  ab  initio  calculations.12  While  MNDO  is  known  to  have 
difficulties  with  compounds  containing  nitro  groups,  ’  the 

M1NDO/3  average  absolute  errors  for  these  compounds  is  essentially 
the  same  as  for  all  CHON  compounds  (approximately  9  kcal/mol).1 
Consequently,  prior  to  the  development  of  AMI,  MIND0/3  was  the 
method  of  choice  in  studies  of  nitro -containing  compounds.1,5 
However,  MNDO  and  MIND0/3  both  predict  the  nitro  group  of  C-l-NP  and 
the  C-l-NP*  to  be  perpendicular  to  the  plane  of  the  alkene.  AMI 
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MINDO/3 . 1,11,13  Since  conjugation  of  the  nitro  group  with  the  alkene 
n  bond,  as  well  as  the  positions  of  the  oxygen  atoms  in  the  plane  of 
the  alkene  as  opposed  to  out  of  the  plane,  could  possibly  have  a 
significant  effect  on  calculated  reaction  profiles,  AMI  was  used  for 
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this  study.  Due  to  the  unusually  large  errors  in  the  AMI  heats  of 
formation  of  NO,  NO^  and  the  larger  than  average  errors  for  N0+, 
N0z+,  and  OH,  experimental  heats  of  formation  for  these  species  are 
used  in  calculating  heats  of  reactions  for  reactions  involving  them. 


All  calculations  were  carried  out  using  AMI  as  incorporated 
in  the  AMPAC  package  of  computer  programs14  modified  to  include 
additional  gradient  norm  minimization  routines.  All  geometries  were 
optimized  by  minimizing  the  energy  with  respect  to  all  geometrical 
variables  using  the  DFP  method.15  Approximate  transition  states 
located  by  reaction  path  calculations  were  optimized  by  minimizing 
the  gradient  norm.18  18  True  minima  and  transition  states  were 
characterized  by  the  presence  of  exactly  zero  or  one  negative 
eigenvalue,  respectively,  in  corresponding  Hessian  matrices.163 


All  calculations  involving  radicals  were  initially  attempted 
in  the  restricted  Hartree-Fock  (RHF)  formalism  with  the 
half -electron  (HE)  approximation. 19  Reaction  path  calculations  for 
many  of  the  reactions  studied  gave  discontinuous  results  using  this 
method.  The  spin-unrestricted  version  of  AMI  was  also  tried,  but 
did  not  eliminate  the  discontinuities.  As  a  result,  2x2 
configuration  interaction  (Cl)  (involving  the  highest  doubly 
occupied  MO  (HDMO)  and  the  singly  occupied  MO  (SOMO))  or  3x3  Cl 
(involving  the  HDMO,  SOMO  and  lowest  unoccupied  MO  (LUMO))  was 
required  to  obtain  continuous  reaction  path  results .  In  all  cases 
but  one,  the  Cl  served  only  to  reach  regions  of  the  geometric 
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surface  that  were  unattainable  using  single  determinantal  wave 
functions.  In  these  cases,  the  approximate  transition  states 
eventually  located  showed  no  mixing  of  configurations ,  and  were 
optimized  satisfactorily  with  the  RHF/HE  method.  The  exception  was 
the  dissociation  of  a  hydroxyl  radical  from  the  product  of  hydrogen 
transfer  to  the  nitro  group,  for  which  3x3  Cl  was  required  to  locate 
and  optimize  the  transition  state.  Results  are  given  in  Table  2. 

Results 

K.  Ionization  of  C-l-NP.  The  heat  of  formation  of  the 
product  of  vertical  ionization  of  C-l-NP  was  calculated  as  256.8 
kcal/mol,  while  that  of  the  optimized  C-l-NP*  was  248.8  kcal/mol. 
The  only  significant  geometrical  change  in  the  relaxation  of 
C-l-NP*  was  the  equalization  of  the  two  carbon-carbon  bonds  (1.34  A 
and  1.47  A  in  C-l-NP,  1.43  A  and  1.42  A  in  C-l-NP*).  This  was 
consistent  with  a  singly  occupied  carbon- carbon  n  orbital  in  the 
cation,  as  the  eigenvectors  indicated. 

B.  Dissociation  of  NO  and  N0+  from  C-l-NP*.  Initial 

2  2 

attempts  to  locate  a  transition  state  (TS)  for  the  dissociation  of 
N02  or  N0z+  using  RHF/HE  calculations  and  treating  the  R-N02  bond  as 
the  reaction  coordinate  resulted  in  discontinuous  results.  Analysis 
of  the  vectors  showed  a  reversal  of  the  HDMO  and  SOMO  at  the  point 
of  discontinuity.  Using  2x2  Cl,  an  approximate  TS  was  located  which 
had  no  mixing  of  configurations.  The  TS  then  optimized  with  the 
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Table  2.  AMI  calculated  heats  of  formation  (AH^),  heats  of  activation  (AH  ) 

and  heats  of  reaction  (AH  )  for  the  mass  spectral  fragmentation  of 

rxn 

cis-l-nitropropene  (kcal/mol). 


AH 

f 

Raaction 

Reactant 

IS 

Product 

ah’* 

AH 

rxn 

C-l-NP  -  C-l-NP* 

5.9 

248. 8“ 

242.9 

C-l-NP*  -  C  H  +  +  NO 

3  5  2 

248.8* 

266.2 

b 

259.1 

17.4 

b 

10.3 

-  C  H  +  NO* 

3  5  2 

b 

287.6 

b 

38.8 

H  C  ONo""l* 

+  3  \  / 

C-1-NP-  -  C-C 

/  \ 

H  H 

248.8* 

268.0 

189.8 

19.2 

>59.8 

3 

3  -  C  H  0  +  N0+ 

3  5 

189.8 

- 

b 

222.0 

- 

b 

32.2 

-  C  H  0*  +  NO 

3  5 

- 

b 

212.2 

- 

b 

22.4 

HO  1+ 

\ 

a  c  no 

+  2  \  / 

C-l-NP-  -  7C-C^ 

248.8* 

239.0 

229.9 

10.2 

-18.9 

H  H 

4 

*  -  C  H  NO*  +  OH 

3  4 

229.9 

255. 3° 

b 

244.9 

c 

25.4 

b 

15.0 

5 

°\  “ 1+ 

H  C  NOH  1  • 

+  3  \  ✓ 

C-l-NP'  -  C-C 

/ 

B 

248.8* 

298.6 

246.5 

49.8 

-2.3 

o  — 1+ 

✓  N  ,0l 

H  C  N 

♦  2  ,  / 

C-l-NP-  -  »  '  +H 

C-C 

248.8* 

286.5 

272.9 

37.7 

24.1 

a  *  b  9  n 

AH  for  the  optimized  C-l-NP- .  The  experimental  AH  values  for  NO  (7.9 

kcal/mol) ,  NO  (231.3  kcal/mol),  NO  (21,6  kcal/mol),  NO  (236.6  kcal/mol)  and 
2  ,  c 

OH  (9.3  kcal/mol)  are  used  in  these  calculations.  Ref.  21. 
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RHF/HE  method  at  1.8  A,  with  a  C-C-C-N  dihedral  angle  of  12°  (as 
opposed  to  0°  in  C-l-NP*).  The  Cl  wavefunction  evidently  enabled 
the  breaking  of  the  planar  symmetry  and,  as  a  result,  a  smooth 
transition  of  orbitals  from  reactants  to  products.  The  lone 
negative  frequency  in  the  force  constant  analysis  was  -823  cm  ^ ,  and 
the  components  of  the  corresponding  normal  vector  primarily  involved 
motion  of  C1  and  the  hydrogens  bonded  to  C1  and  C2 .  This  motion 
would  ultimately  restore  planarity  to  the  hydrocarbon  fragment  and 
increase  the  bond  angle  of  the  hydrogen  bonded  to  C1  from  123°  to 
near  180°  in  the  C^H^  product.  The  enthalpy  of  activation  and 
enthalpy  of  reaction  for  the  formation  of  C  H  +  +  NO  were  17.4  and 

3  5  2 

10.3  kcal/mol,  respectively.  Dissociation  to  produce  and  N02+ 
is  predicted  to  be  30  kcal/mol  more  endothermic . 

C.  Nitro  to  nitrite  rearrangement.  The  approximate  TS  was 
located  with  RHF/HE  calculations  by  treating  the  C^N-O  bond  angle 
for  the  oxygen  closest  to  the  methyl  group  as  the  reaction 
coordinate.  In  the  optimized  TS,  Figure  2,  the  nitro  group  is 
virtually  perpendicular  to  the  plane  of  the  alkene.  The  incipient 
C-0  bond  is  1.7  A,  while  the  corresponding  N-0  bond  has  stretched 
only  from  1.20  A  in  C-l-NP-  to  1.26  A.  In  the  product,  3,  the  C-0 
bond  is  1.30  A  and  the  0-N0  bond  is  extremely  long,  1.69  A.  The 
negative  frequency  of  the  force  matrix  is  -705  cm  1 ,  the  components 
of  the  vector  primarily  involving  motion  of  Cl  and  the  N02  group, 
although  there  is  also  significant  participation  of  the  hydrogen 
bonded  to  C2.  The  rearrangement  is  predicted  to  be  very  exothermic 
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Figura  2.  Transition  stats  for  ths  nitro  to  nitrits  rsarrangeasnt  cf  C-l-NP- .  The  0  -N-C  -C 
o 

dihsdrsl  angla  is  94.4  .  In  this  and  ths  following  figurss,  bond  lsngths  ars  in  angstroms  and 
bond  anglss  ars  in  dsgrsss.  Atoms  not  Labslsd  ars  hydrogsns. 

(AH  -  -59.8  kcal/mol)  and  the  enthalpy  of  activation,  19.2 

rxn 

kcal/mol,  is  similar  to  that  for  N0z  dissociation.  The  strain 

associated  with  the  three  membered  cyclic  TS  for  the  rearrangement 

should  decrease  the  probability  of  the  rearrangement  occurring 

relative  to  loss  of  NO  .  When  it  does  occur,  however,  the  80 

2 

kcal/mol  of  reverse  activation  energy  stored  in  the  product  makes 
rapid  decomposition  of  3  likely. 

D.  Dissociation  of  NO  and  N0+  from  3.  No  TS  could  be 
located  for  dissociation  of  NO  or  N0+  from  3.  This  is  reasonable 
considering  there  is  very  little  geometry  change  expected  in  the 
fragments  due  to  the  dissociation.  The  22  and  32  kcal/mol  required 
to  dissociate  NO  and  N0+,  respectively,  are  certainly  insufficient 


to  prevent  the  fragmentations  given  the  vibrational  energy  stored  in 
3  as  a  result  of  the  nitro  to  nitrite  rearrangement.  At  first 
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Figure  3.  Transition  state  for  hydrogen  transfer  from  C3  of  C-l-NP-  to  01  .  The  H-C3-C2-C1 
o  112  o 

dihedral  angle  is  -30. *  ,  and  the  0  -N-C  -C  dihedral  angle  is  8.6  . 
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glance  it  seems  surprising  the  fragmentation  to  C^H^O  +  NO  requires 
less  energy  than  the  fragmentation  to  C^O  +  N0+,  since  there  is 
only  a  very  small  m/z  57  peak  (relative  abundance  3.2)  while  there 
is  a  sizeable  m/z  30  peak  (relative  abundance  17.1).  This  will  be 
discussed  further  in  the  Discussion  Section. 


E.  Hydrogen  transfer  from  the  methyl  group  to  the  nitro 
group.  As  with  the  N0z  dissociation  from  C-l-NP*,  2x2  Cl 
calculations  were  required  to  eliminate  discontinuities  observed  for 
RHF/HE  calculations  treating  the  H-0  distance  as  the  reaction 
coordinate.  The  approximate  TS  located  was  then  optimized  with  the 
RHF/HE  method,  as  there  was  no  mixing  of  configurations  at  the 
predicted  TS.  The  H-0  distance  in  the  optimized  TS,  Figure  3,  is 
1.62  A,  while  the  H-C  bond  length  increased  from  1.14  to  1.24  A  and 
the  N-0  bond  length  from  1.20  to  1.25  A.  The  oxygen  atom  and 


hydrogen  atom  were  0.5  and  0.6  A  out  of  the  C-C-C  plane, 

respectively.  The  negative  frequency  of  the  force  matrix  was  -842 

cm  \  the  motion  of  the  normal  vector  being  composed  of  motion  of 

the  transferring  hydrogen,  as  well  as  adjustment  of  the  other 

hydrogens  bonded  to  C3.  The  10.2  kcal/mol  activation  barrier  is 

lower  than  those  for  NO  dissociation  and  nitro  to  nitrite 

2 

rearrangement.  The  formation  of  the  TS  does  limit  the  motion  of  the 
methyl  rotor,  resulting  in  an  unfavorable  entropic  contribution  to 
the  formation  of  the  TS.  At  least  30  kcal/mol  of  reverse  activation 
energy  is  stored  in  the  product,  4,  as  a  result  of  the  hydrogen 
transfer. 

F.  Dissociation  of  OH  from  4.  A  3x3  Cl  calculation  was 

«»• 

required  to  locate  and  optimize  the  TS  for  the  dissociation  of  the 
OH  radical  from  4.  In  the  TS ,  the  N-OH  bond  length  was  2.0A  and 
the  0-N-C-C  dihedral  was  6°,  as  compared  to  1.3  A  and  1° , 
respectively,  in  4.  The  negative  frequency  in  the  force  matrix  was 
-443  cm  \  and  the  normal  vector  was  composed  primarily  of  the 
motion  of  the  OH  group.  The  25.4  kcal/mol  activation  barrier21  was 
more  than  compensated  for  by  the  reverse  activation  energy  from  the 
rearrangement,  which  was  stored  as  vibrational  energy  in  4. 

G.  Hydrogen  transfer  from  C1  to  the  nitro  group.  While  the 
presence  of  the  m/z  70  peak  in  the  C-l-NP  spectrum  is  interpreted  as 
a  sign  of  interaction  between  the  methyl  and  nitro  groups  during 
fragmentation,  these  same  peaks  could  result  from  hydrogen  transfer 
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Figure  4.  Transition  stata  for  hydrogen  transfer  from  C  of  C-l-NP*  to 
o  12  3  o 

dihedral  angle  la  132.7  ,  while  the  N-C  -C  -C  dihedral  angle  ia  -40.3  . 
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from  C1  to  the  nitro  group.  This  transfer  could  occur  in  both 
C-l-NP  or  T-l-NP.  If  the  enthalpic  requirement  for  this  reaction 
was  similar  to  that  for  hydrogen  transfer  from  the  methyl  group,  a 
different  explanation  might  be  required  for  the  differences  in  the 
two  spectra.  As  with  the  previous  hydrogen  transfer,  2x2  Cl  was 
required  to  locate  the  TS ,  but  it  was  optimized  satisfactory  with 
the  RHF/HE  method.  The  TS ,  Figure  4,  was  extremely  strained.  The 
negative  frequency  was  -2244  cm  ^  and  corresponded  entirely  to  the 
motion  of  the  migrating  hydrogen.  The  high  activation  barrier,  49.8 
kcal/mol,  plus  the  strain  in  the  transition  state  make  this 
rearrangement  highly  improbable,  supporting  the  previous 
interpretation  of  the  m/z  70  peak. 


Figur#  5.  Transition  state  for  the  loss  of  a  hydrogen  atom  from  C  of  C-l-NP- 
13  2  1  o 

0  -C  -C  -C  dihedral  angle  is  0.0  . 


The 


H.  Dissociation  of  a  hydrogen  atom  from  the  methyl  carbon. 
It  is  interesting  to  note  the  significant  m/z  86  peak  (relative 
abundance  11.3)  in  the  C-l-NP  spectrum  and  the  negligible  peak 
(relative  abundance  1.9)  in  the  T-l-NP  spectrum.  When  the  C3-0 
distance  was  treated  as  the  reaction  coordinate  and  incrementally 
decreased,  a  hydrogen  atom  eventually  dissociated  from  the  methyl 
group  and  a  cationic  nitrogen  analog  to  a  lactone  was  formed.  The 
same  product  resulted  when  the  length  of  a  C-1-NP+  H-C3  bond  was 
incrementally  increased.  The  TS  for  this  reaction  was  located  by 
starting  with  the  product  of  the  dissociation  and  incrementally 
decreasing  the  H---C3  distance.  The  RHF/HE  method  was  used  to 
locate  and  optimize  the  TS,  Figure  5.  The  H-C3  distance  in  the  TS 
was  1.54  A,  while  the  incipient  C3-0  bond  was  1.63  A.  The  37.7 


kcal/mol  activation  is  significantly  larger  than  :hose  of  the  other 
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occurring  reactions,  and  formation  of  this  TS  again  restricts  the 
motion  of  the  methyl  rotor. 


Discussion 

Electron  impact  ionization  produces  molecular  ions  with 
internal  energies  ranging  from  zero  to  over  20  eV,  the  average  being 
a  few  eV.22  According  to  the  quasi -equilibrium  theory,23  the 
statistical  theory  generally  accepted  as  describing  the 
fragmentation  in  mass  spectrometers,22,24  this  internal  energy  is 
rapidly  and  randomly  distributed  among  the  internal  degrees  of 
freedom  in  the  ground  state  of  the  molecular  ion.  Fragmentation  is 
then  statistical,  specific  fragmentations  occurring  only  when  the 
nuclei  are  in  the  proper  configuration  and  sufficient  vibrational 
energy  is  concentrated  in  the  necessary  degrees  of  freedom.23 
Arguments  involving  fragmentation  from  isolated  electronic  states 
and  nonrandom  distribution  of  vibrational  energy  are  sometimes  used 
to  explain  apparently  nonstatistical  fragmentation. 8,25  Nitromethane 
and  nitrobenzene  are  two  compounds  for  which  these  arguments  have 
been  proposed.*  It  is  difficult  to  eliminate  alternative 
explanations  such  as  noninterconverting  structural  isomers,8'26 
however,  and  incomplete  randomization  has  not  been  proven  directly.8 
Consequently,  it  is  assumed  C-1-NP+  fragments  statistically  from 
its  ground  electronic  state.  Additionally,  while  there  is 
experimental  evidence  that  some  fragmentations  formally  described  as 
direct  bond  cleavages  actually  involve  more  extensive  intramolecular 
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interactions,  including  dissociation  of  N02  from  nitrobenzene, 
all  dissociations  in  this  study  are  assumed  to  involve  loose 
transition  states. 

Attempts  to  rationalize  low  m/z  peaks  are  often  futile 
because  they  are  likely  to  correspond  to  products  of  secondary 
fragmentation,  and  are  also  more  likely  to  involve  extensive 
rearrangements.  The  exception  is  when  there  is  an  obvious, 
straight-forward  path  to  the  low  m/z  fragment,  such  as  m/z  30  (C^O 
+  NO+)  and  m/z  41  (C  H  +  +  NO  ) ,  and  even  then  there  is  no  guarantee 
this  path  accounts  for  the  full  peak  abundance.  Higher  m/z 
fragments,  while  still  possibly  involving  rearrangement,  will  more 
probably  correspond  to  products  of  fragmentation  of  the  molecular 
ion  or  its  isomer.  The  m/z  86  and  m/z  70  peaks,  which  are  present 
in  the  C-l-NP  spectrum  and  virtually  absent  in  the  T-l-NP  spectrum, 
fall  in  this  latter  category,  and  are  easily  rationalized  from  AMI 
results  on  the  basis  of  direct  participation  of  the  methyl  group  in 
the  fragmentation  of  C-l-NP*.  The  AMI  results  also  explain  certain 
similarities  in  the  two  spectra. 

In  both  spectra,  there  are  large  m/z  41  peaks  corresponding 
to  fragmentation  to  C3H  +  +  N0z>  and  small  m/z  46  peaks, 

corresponding  to  fragmentation  to  +  N02+.  This  would  certainly 
be  expected  on  the  basis  of  AMI  calculations,  which  indicate  it  is 
more  favorable  for  the  hydrocarbon  fragment  to  carry  the  charge  by 
30  kcal/mol.  While  the  formation  of  CH+  +  NO  from  the 
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nitrobenzene  molecular  ion  apparently  occurs  with  no  reverse 
activation  energy,4  there  is  a  small  reverse  activation  energy  for 
the  C-l-NP*  fragmentation.  This  is  reasonable  since  there  is 
significant  geometric  relaxation  of  the  C3H5+  fragment  due  to 

fragmentation,  whereas  small  changes  would  be  expected  in  the 

relatively  rigid  phenyl  cation. 

On  the  basis  of  the  calculated  enthalpy  of  activation  for 
the  nitro  to  nitrite  rearrangement,  significant  rearrangement  should 
occur,  and,  as  indicated  by  the  size  of  the  m/z  30  peak,  does.  The 
negative  frequency  of  the  TS  is  relatively  low,  indicating  the 
strain  in  the  three  membered  cyclic  TS  is  not  excessively  large  and 
is  insufficient  to  make  the  rearrangement  uncompetitive. 

The  AMI  results  for  fragmentation  of  3  to  +  NO  and 

C3H50  +  N0+  do  seemingly  conflict  with  the  presence  of  m/z  30  peaks 

and  the  absence  of  m/z  57  peaks  in  the  C-l-NP  and  T-l-NP  spectra. 

These  results  are,  however,  perfectly  consistent.  The  nitro  to 
nitrite  rearrangment  occurs  with  a  large  reverse  activation  energy, 
resulting  in  approximately  80  kcal/mol  of  vibrational  energy  being 
stored  in  3.  Upon  fragmention  of  3,  the  excess  energy  is 

partitioned  between  the  two  fragments ,  with  the  larger  fragment 
expected  to  carry  off  more  of  the  energy  (degree  of  freedom 
effect23).  As  a  result,  the  C3Hs0+  ion  will  likely  fragment 
further.  One  possibility  is  that  it  loses  a  water  molecule,  thereby 
contributing  to  the  large  m/z  39  peaks  in  both  spectra.  On  the 


being  stored  in  4.  This  is  more  than  enough  energy  to  result  in 

21 

subsequent  decomposition  to  5  +  OH,  and  the  unfavorable  entropy 
factor  associated  with  the  restriction  of  the  methyl  rotor  in  the  TS 
does  not  outweigh  the  relatively  favorable  enthalpic  factor.  These 
results,  plus  the  high  energy  and  tight  TS  required  for  hydrogen 
transfer  from  C1  to  the  nitro  group,  clearly  explain  the  presence  of 
the  m/z  70  peak  in  the  C-l-NP  spectrum  and  its  absence  in  the  T-l-NP 
spectrum. 

Given  the  energetic  requirements  for  the  hydrogen  transfer 
and  hydroxyl  radical  dissociation  versus  those  for  the  dissociation 
of  NO^,  it  is  even  a  little  surprising  the  m/z  70  peak  is  less  than 
half  the  size  of  the  m/z  41  peak  in  the  C-l-NP  spectrum.  One  reason 
may  be  the  restriction  of  the  methyl  rotor  in  the  TS  for  hydrogen 
transfer.  Additionally,  there  is  another  route  to  the  production  of 
an  m/z  41  peak,  namely  loss  of  N0z  from  3.  The  heat  of  formation 
for  the  TS  of  this  process  would  be  expected  to  be  similar  to  that 
for  the  loss  of  N0z  from  C-l-NP*.  This  is  probably  a  minor 
process,  however,  since  loss  of  NO  and  N0+  is  much  more  heavily 
favored  energetically,  as  would  be  anticipated  from  the  long  0-N0 


bond.  In  any  case,  It  is  unlikely  these  possibilities  provide  the 
full  explanation. 


On  the  other  hand,  a  major  difference  between  the  C-l-NP  and 
T-l-NP  spectra  still  unaccounted  for  is  the  m/z  42  peak.  Its  large 
size  in  the  C-l-NP  spectrum  and  virtual  absence  in  the  T-l-NP 
spectrum  indicates  its  generation  in  some  way  involves  interaction 
between  the  methyl  and  nitro  groups.  Only  three  combinations  of 
atoms  present  in  C-l-NP  add  up  to  an  m/z  42:  C^O*.  C^N*  and 
CNO+.  A  reasonable  source  of  the  m/z  42  ion  consistent  with  the 
above  data,  then,  would  be  fragmentation  of  5  to  ethylene  and  CNO+. 
While  the  total  heat  of  formation  for  these  products  is  predicted  to 
be  quite  high  <CNO+  342.4  kcal/mol  (AMI),  CH  12.5kcal/mol 
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(exptl20) ,  it  is  still  reachable,  considering  the  internal  energy 
deposited  in  the  molecular  ion  during  ionization.  Other 
explanations  for  the  m/z  42  peak  are  either  inconsistent  with  the 
contrast  between  the  C-l-NP  and  T-l-NP  spectra  or  involve 
significantly  more  rearrangement. 

Finally,  the  dissociation  of  a  hydrogen  atom  from  the  methyl 
group  with  simultaneous  formation  of  the  five  membered  cyclic  ion  is 
also  consistent  with  the  two  spectra,  i.e  the  much  larger  m/z  86 
peak  in  the  C-l-NP  spectrum  as  opposed  to  the  T-l-NP  spectrum. 
While  the  enthalpy  of  activation  and  endothermicity  for  this 
reaction  is  significantly  larger  than  the  other  occurring 
intramolecular  processes,  they  are  less  than  would  be  expected  if  a 


simple  C-H  bond  scission  was  occurring  without  the  assistance  of 
other  intramolecular  interactions . 

Conclusions 

Differences  in  the  mass  spectra  of  C-l-NP  and  T-l-NP  imply 
there  is  significant  involvement  of  the  methyl  group  in  the 
fragmentation  of  C-l-NP*.  The  results  of  AMI  calculations  for 
several  possible  fragmentation  paths  support  this  conclusion  and 
appear  to  provide  a  good  rationalization  for  the  major  features  of 
the  C-l-NP  spectrum.  Further  experimental  work,  involving  softer 
ionization  and  collision  induced  dissociation,  is  necessary  to 
confirm  the  proposed  looseness  of  the  transition  states  for 
dissociations  and  secondary  fragmentation  pathways. 

Based  on  these  results ,  it  can  be  concluded  that  the  MINDO/3 
and  MNDO  results  of  Turner  and  Davis1  for  possible  thermal 
decomposition  mechanisms  for  C-l-NP  are  probably  also  correct.  As 
noted  by  them,  however,  using  C-l-NP  as  a  model  for  TNT  has  a 
shortcoming  in  that  the  effects  on  the  aromaticity  of  TNT  is  not 
represented  in  calculations  involving  C-l-NP.  This  is  the  most 
likely  reason  the  R-N02  bond  scission,  predicted  to  be  unfavorable 
for  C-l-NP,  was  found  to  be  the  dominant  process  in  the  gas  phase 
pyrolysis  of  substituted  nitrobenzenes . 
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ARCHIVE  OF  SELECTED  GEOMETRIES  FROM  CHAPTER  4 


METHYL  MALONATE,  HYDROGEN  BOND  ON  CARBONYL  FROM  CH3NH3+ 
REACTANT  FOR  DECARBOXYLATION,  -134.3  KCAL/MOL  AMI 
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METHYL  MALONATE, 
TRANSITION  STATE 


HYDROGEN  BOND  ON  CARBONYL  FROM  MENH3+ 
FOR  DECARBOXYLATION,  -123.6  KCAL/MOL  AMI 
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METHYL  MALONATE 
HYDRATED  WITH  5 


(MALONYL-ACP  MODEL) ,  CARBONYL  OXYGEN 
WATER  MOLECULES,  -470.1  KCAL/MOL  AMI 
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METHYL  MALONATE,  HYDROGEN  BOND  ON  CARBONYL 
REACTANT  FOR  DECARBOXYLATION,  -311.1  KCAL/MOL  AMI 
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METHYL  THIOACETATE  +  ENTHIOLATE 

TRANSITION  STATE  OF  CONDENSATION,  -90.2  KCAL/MOL  MNDO 
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METHYL  THIOLACETATE  +  ENTHIOLATE 

PRODUCT  OF  CONDENSATION,  -108.1  KCAL/MOL  MNDO 
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METHYL  ACETATE  +  ENOLATE 

CHARGE-DIPOLE  COMPLEX,  -198.0  KCAL/MOL  AMI 
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METHYL  ACETATE  +  ENOLATE 

TRANSITION  STATE  FOR  CONDENSATION  FROM  CHARGE-DIPOLE  COMPLEX, 
-192.2  KCAL/MOL  AMI 
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METHYL  ACETATE  +  ENOLATE 
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PRODUCT  OF  CONDENSATION,  -199. 

9  KCAL/MOL  AMI 
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1 
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C 
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1 

6 
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8 
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H 
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1 
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12 

6 

1 

C 

1.479715 

1 
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1 
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1 

12 

13 
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15 

12 
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15 

12 

14 

C 

1.419069 

1 

116.933165 

1 

-174.114868 

1 

17 

15 

12 

H 

1.116811 

1 

110.424977 

1 

55.476836 

1 
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17 
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MALONATE  ESTER 

4-  METHYL  AMMONIUM 

ION  +  METHYL 

ACETATE 

COMPLEX,  -240.4 

KCAL/MOL  AMI 

XX 

0 . 000000 

0 

0.000000 

0 

0 . 000000 

0 

0 

0 

0 

0 

1.043215 

1 

0.000000 

0 

0 . 000000 

0 

1 

0 

0 

C 

1.425681 

1 

88.079845 

1 
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0 

2 

1 

0 

H 

1.116864 

1 

110.292107 

1 

-60.853607 

1 

3 

2 

1 

H 

1.118876 

1 

110.103310 

1 

61.997164 

1 

3 

2 

1 

H 

1.119381 

1 

104.059443 

1 

-179.680706 

1 

3 

2 
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1 

8 
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1 

8 
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H 
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1 
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1 

-51.925054 

1 

8 
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1 
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1 

7 

9 

2 

H 

1.115842 

1 
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13 

7 
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13 

7 
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1 
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1 

16 

13 

7 

0 
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1 
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1 
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1 

16 

13 
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C 
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1 
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1 

9.524380 

1 

18 

16 

17 

H 

1.116223 

1 

110.273844 

1 
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1 

19 

18 

16 

H 

1.116288 

1 
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1 

-68.384845 

1 

19 

18 

16 

H 

1.117939 

1 
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1 
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1 

19 

18 

16 

C 
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1 
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0 
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0 

13 

7 
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0 
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23 

13 

16 
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23 
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7 

2 

N 

1.031582 

1 

139.968946 

1 

37.108715 

1 

26 

9 

7 

H 

1.045875 

1 

100.627324 

1 

105.714146 

1 

27 

9 

7 

H 

1.026423 

1 

107.967539 

1 

-116.035089 

1 

27 

26 

28 

C 

1.467055 

1 

111.035620 

1 

122.361936 

1 

27 

28 

29 

H 

1.121976 

1 

109.933648 

1 

-176.579478 

1 

30 

27 

28 

H 

1.123291 

1 

109.457170 

1 

-56.237352 

1 

30 

27 

28 

H 

1.123361 

1 

109.147937 

1 

63.348500 

1 

30 

27 

28 

0 

0 . 000000 

0 

0 . 000000 

0 

0 . 000000 

0 

0 

0 

0 

. 


Page  159 


MALONATE  ESTER  +  METHYL  AMMONIUM  ION  +  METHYL  ACETATE 
TRANSITION  STATE  FOR  DECARBOXYLATION,  -235.2  KCAL/MOL  AMI 


XX 

0.000000 

0 

0.000000 

0 

0 . 000000 

0 

0 

0 

0 

0 
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0 
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0 

1 
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3 

2 
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H 
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1 
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13 
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1 
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0 

ENTHIOLATE  +  METHYL  AMMONIUM  ION  +  METHYL  ACETATE 
PRODUCT  OF  DECARBOXYLATION,  REACTANT  FOR  CONDENSATION, 

-180.1  KCAL/MOL  AMI 
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ENTHIOLATE  +  METHYL  AMMONIUM  ION  +  METHYL  ACETATE 
CONDENSATION  PRODUCT  WITH  DOUBLE  H+  TRANSFER, 


-206. 

.2  KCAL/MOL  AMI 
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1 
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1 
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1 

13 

7 
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18 

16 
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1 
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1 
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1 

19 

18 

16 
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0 
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0 

0.000000 

0 

13 

7 

8 
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0 

114.785398 

0 

77.815034 

0 

23 

13 

16 
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0 

129.404298 

0 

-177.018699 

0 

23 

24 

13 
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1 

108.154233 

1 

43.214092 

1 

9 

7 

2 

2.806837 

1 

150.822977 

1 

-56.984304 

1 

26 

9 

7 

1.002480 

1 

70.370855 

1 

115.097967 

1 

27 

9 

7 

1.002878 

1 

52.102014 

1 

-147.397303 

1 

27 

26 

28 
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1 
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1 
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1 

27 

28 

29 
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1 
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1 

30 

27 

28 
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30 

27 

28 
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27 

28 
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APPENDIX  2 


ARCHIVE  OF  TRANSITION  STATE  GEOMETRIES  FROM  CHAPTER  5 


CIS - 1 -NITROPROPENE  CATION 

RHF/HE  TRANSITION  STATE  FOR  N02  LOSS,  266.2  KCAL/MOL  AMI 


c 

0.000000 

0 

0.000000 

0 

0.000000 

0 

0 

0 

0 

c 

1.472561 

1 

0.000000 

0 

0 . 000000 

0 

1 

0 

0 

C 

1.310066 

1 

130.892582 

1 

0.000000 

0 

2 

1 

0 

N 

1.816906 

1 

115.211682 

1 

11.844782 

1 

3 

2 

1 

0 

1.158931 

1 

115.816566 

1 

78.199726 

1 

4 

3 

2 

0 

1.168146 

1 

107.309388 

1 

-103.087237 

1 

4 

3 

2 

H 

1.131380 

1 

113.603362 

1 

-163.782231 

1 

2 

3 

4 

H 

1.094209 

1 

145.231752 

1 

-169.728677 

1 

3 

2 

4 

H 

1.123058 

1 

110.483069 

1 

64.775317 

1 

1 

2 

3 

H 

1.121958 

1 

111.050414 

1 

-175.502058 

1 

1 

2 

3 

H 

1.122378 

1 

110.366934 

1 

-55.665578 

1 

1 

2 

3 

0 

0.000000 

0 

0.000000 

0 

0 . 000000 

0 

0 

0 

0 

CIS -1 -NITROPROPENE  CATION 

RHF/HE  TRANSITION  STATE  FOR  THE  NITRO  TO  NITRITE 
REARRANGEMENT,  268.0  KCAL/MOL  AMI 


C 

0.000000 

0 

0.000000 

0 

0.000000 

0 

0 

0 

0 

c 

1.446441 

1 

0.000000 

0 

0.000000 

0 

1 

0 

0 

c 

1.411156 

1 

126.767102 

1 

0.000000 

0 

2 

1 

0 

N 

1.478678 

1 

123.397025 

1 

-15.173841 

1 

3 

2 

1 

0 

1.151599 

1 

153.169802 

1 

-87.847290 

1 

4 

3 

2 

0 

1.265076 

1 

76.137635 

1 

94.360088 

1 

4 

3 

2 

H 

1.113172 

1 

115.064054 

1 

165.623437 

1 

2 

3 

4 

H 

1.118437 

1 

123.083732 

1 

-169.411905 

1 

3 

2 

4 

H 

1.130221 

1 

109.803527 

1 

119.694335 

1 

1 

2 

3 

H 

1.131032 

1 

109.266842 

1 

-125.002679 

1 

1 

2 

3 

H 

1.117745 

1 

114.834557 

1 

-3.245965 

1 

1 

2 

3 

0 

0.000000 

0 

0.000000 

0 

0.000000 

0 

0 

0 

0 
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CIS-l-NITROPROPENE  CATION 

RHF/HE  TRANSITION  STATE  FOR  H  TRANSFER  FORM  C3  TO  0, 


259 

.0  KCAL/MOL 

AMI 

C 

0.000000 

0 

0.000000 

0 

0.000000 

0 

0 

0 

0 

C 

1.447477 

1 

0 . 000000 

0 

0.000000 

0 

1 

0 

0 

C 

1.349242 

1 

123.593025 

1 

0 . 000000 

0 

2 

1 

0 

N 

1.482576 

1 

119.901972 

1 

7.946509 

1 

3 

2 

1 

0 

1.180147 

1 

121.468956 

1 

-172.173900 

1 

4 

3 

2 

0 

1.248468 

1 

118.479657 

1 

8.574945 

1 

4 

3 

2 

H 

1.117041 

1 

118.908619 

1 

-171.004935 

1 

2 

3 

4 

H 

1.109687 

1 

127.434873 

1 

-179.592901 

1 

3 

2 

4 

H 

1.124623 

1 

114.874983 

1 

78.477856 

1 

1 

2 

3 

H 

1.119693 

1 

115.698610 

1 

-146.509560 

1 

1 

2 

3 

H 

1.617934 

1 

114.492292 

1 

-1.953327 

1 

6 

4 

3 

0 

0.000000 

0 

0 . 000000 

0 

0 . 000000 

0 

0 

0 

0 

PRODUCT  OF  H  TRANSFER  FROM  C3 

OF 

C-1-NP+  TO  0 

3X3 

Cl  TRANSITION 

STATE  FOR  LOSS 

OF  OH,  255.3 

KCAL/MOL  AMI 

•  C 

0.000000 

0 

0.000000 

0 

0 . 000000 

0 

0 

0 

0 

C 

1.373191 

1 

0.000000 

0 

0.000000 

0 

1 

0 

0 

C 

1.400414 

1 

128.829508 

1 

0 . 000000 

0 

2 

1 

0 

N 

1.361015 

1 

127.612797 

1 

5.402652 

1 

3 

2 

1 

0 

1.137269 

1 

139.111242 

1 

-178.165153 

1 

4 

3 

2 

0 

2.019259 

1 

116.489436 

1 

5.759621 

1 

4 

3 

2 

H 

1.114726 

1 

112.795416 

1 

-175.352421 

1 

2 

3 

4 

H 

1.120051 

1 

118.841171 

1 

179.556314 

1 

3 

2 

4 

H 

1.102970 

1 

124.627890 

1 

2.756380 

1 

1 

2 

3 

H 

1.106078 

1 

119.822026 

1 

-178.239945 

1 

1 

2 

3 

H 

0.954933 

1 

109.620139 

1 

179.395724 

1 

6 

4 

3 

0 

0.000000 

0 

0.000000 

0 

0 . 000000 

0 

0 

0 

0 

CIS • 1 -NITROPROPENE  CATION 
RHF/HE  TRANSITION  STATE  FOR  H 
298 . 6  KCAL/MOL  AMI 


0.000000 

1.456441 

1.356647 

1.451867 

1.250787 

1.172156 

1.119283 

1.434974 

1.132306 

1.127726 

1.119358 

0.000000 


CIS - 1 -NITROPROPENE 
RHF/HE  TRANSITION 
286.5  KCAL/MOL  AMI 


STATE  FOR  H  TRANSFER  FROM  Cl  TO  0, 

0.000000  0  0.000000  0  0  0  0 

0.000000  0  0.000000  0  100 

124.908753  1  0.000000  0  210 

140.806359  1  -40.275652  1  321 

103.501132  1  177.694610  1432 

131.397423  1  -3.774432  1432 

116.889250  1  147.516712  1  234 

85.698827  1  -2.154217  1  543 

107.751755  1  124.966864  1  123 

110.622908  1  -119.132544  1  123 

114.209343  1  4.240303  1123 

0.000000  0  0.000000  0  0  0  0 


CATION 

STATE  FOR  LOSS  OF  H  FROM  C3, 


c 

0.000000 

0 

0.000000 

0 

0.000000 

0 

0 

0 

0 

c 

1.499151 

1 

0.000000 

0 

0.000000 

0 

1 

0 

0 

c 

1.353726 

1 

111.995843 

1 

0.000000 

0 

2 

1 

0 

N 

1.480957 

1 

108.832258 

1 

-0.009507 

1 

3 

2 

1 

0 

1.166514 

1 

129.384838 

1 

179.987263 

1 

4 

3 

2 

0 

1.633695 

1 

99.008836 

1 

0.025864 

1 

1 

2 

3 

H 

1.101695 

1 

125.762584 

1 

180.012782 

1 

2 

3 

4 

H 

1.101610 

1 

133.944226 

1 

-179.981520 

1 

3 

2 

4 

H 

1.129910 

1 

117.671652 

1 

101.318940 

1 

1 

2 

3 

H 

1.537874 

1 

97.108168 

1 

179.917626 

1 

1 

2 

3 

H 

1.129771 

1 

117.649361 

1 

-101.301370 

1 

1 

2 

3 

0 

0.000000 

0 

0.000000 

0 

0.000000 

0 

0 

0 

0 

